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ABSTRACT 

The structure of relativistic radiation mediated shocks (RRMS) propagating into a cold electron-proton 
plasma is calculated and analyzed. A qualitative discussion of the physics of relativistic and non relativis- 
tic shocks, including order of magnitude estimates for the relevant temperature and length scales, is presented. 
Detailed numerical solutions are derived for shock Lorentz factors T u in the range 6 < T„ < 30, using a novel 
iteration technique solving the hydrodynamics and radiation transport equations (the protons, electrons and 
positrons are argued to be coupled by collective plasma processes and are treated as a fluid). The shock tran- 
sition (deceleration) region, where the Lorentz factor V drops from T„ to ~ 1 , is characterized by high plasma 
temperatures T ~ Tm e c 2 and highly anisotropic radiation, with characteristic shock-frame energy of upstream 
and downstream going photons of a few x , 



and ~ T 2 m e c 2 , respectively. Photon scattering is dominated 
by e ± pairs, with pair to proton density ratio reaching « 10 2 r„. The width of the deceleration region, in terms 
of Thomson optical depths for upstream going photons, is large, At ~ T 2 (At ~ 1 neglecting the contribution 
of pairs) due to Klein Nishina suppression of the scattering cross section. A high energy photon component, 
narrowly beamed in the downstream direction, with a nearly flat power-law like spectrum, vl v ex ir, and an 
energy cutoff at ~ T 2 m e c 2 carries a fair fraction of the energy flux at the end of the deceleration region. An 
approximate analytic model of RRMS, reproducing the main features of the numerical results, is provided. 
Subject headings: shock waves — radiation mechanisms: nonthermal — gamma-rays: bursts 



1. INTRODUCTION 

Radiation mediated shocks (RMSs) are shocks in which the 
downstream (DS) energy density is dominated by radiation 
rather than by particle thermal energy, and in which the fast 
upstream (US) plasma approaching the shock is decelerated 
by scattering of photons, generated in the DS and propagating 
into the US, by the fast US electrons. RMS are expected to oc- 
cur in a variety of astrophysical flows. The shock waves prop- 
agating through, and expelling, the envelopes of massive stars 
undergoing core collapse supernova explosions, are non rela- 
tivistic (NR) RMS ( Weavedl 19761) . Relativistic RMS (RRMS) 
may play an important role in, e.g., gamma-ray bursts, trans- 
relativistic suprenovae, and pulsar accretion flows. 
[1] Gamma Ray Bursts (GRBs). Within the fram ework of 
the collapsar model of GRBs (e.g. IWooslevlll993l) . a highly 
relativistic jet driven by the collapsed core of a massive star 
penetrates through the stellar envelope. The shock that de- 
celerates the jet is expected to be a highly relativistic RMS 
(Mes zaros & Waxmanl200UlAloy et alJl2000h . 
[2] Trans Relativistic SNe. Several recent SN events, that were 
identified in very early stages of the explosion, have been 
shown to deposit a significant fraction, ~ 1%, of the explosion 
energy in mildly relativistic, 7/3 > 1, ejecta (Soderberg et al. 
2006, and references therein). The existence of mildly rela- 
tivistic ejecta components suggests that a mildly relativistic 
RMS shock traversed the outer envelope of the progenitor. 
[3] Pulsar accretion flows. Accretion onto the polar cap of 
a pulsar is expected to produce a mildly relativistic RMS 
which is appro ximately stationary in the neutron star frame 
(Burnard et al. 1991; Becker 1988, and references therein). 

NR RMS were studied in detail in dWeaverlll976 ). describ- 
ing photon propagation using the diffusion approximation and 
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describing the radiation field using two parameters, photon ef- 
fective temperature and density. These approximations hold 
for slow shocks, v/c < 0.2, for which relativistic effects are 
negligible and the Thomson optical depth of the shock de- 
celeration region is large, ~ c/v (ensuring that the radiation 
field is nearly isotropic and that the photons are in Comp- 
ton equilibrium). The NR approximations do not hold for 
faster shocks. For such shocks, relativistic effects (such as 
pair production and relativistic corrections to the cross sec- 
tions of radiative processes) are important, and the radiation 
field becomes highly anisotropic. 

A simplified solution for the structure of RRMS, neglect- 
ing pair production, photon production and relativistic cor- 
rections, was derived by Levinson & Bromberg (2008). This 
solution may be applicable only in cases where the US plasma 
holds a significant photon density, which keeps the plasma 
at low temperatures throughout the shock, much lower than 
those obtained in a self-consistent solution where the pho- 
ton density vani shes at US infinity. In a preceding paper 
dKatzet al.l 2010) we derived a simple approximate analytic 
model for the structure of radiation mediated shoc ks. This 
model accurately reproduces the numerical results of Weaver 
d 19761) for v/c < 0.2, and provides an approximate descrip- 
tion of the shock structure at larger velocities, v/c — > 1 . We 
confirmed that at shock velocities v/c > 0.1 the shock tran- 
sition region is far from thermal equilibrium, with electrons 
and photons (and positrons) in Compton (pair) equilibrium at 
temperatures T s significantly exceeding the far downstream 
temperature. We have found that T s > lOkeV is reached 
at shock velocities v/c « 0.2, and that at higher velocities, 
v/c > 0.6, the plasma is dominated in the transition region 
by e ± pairs and 60 keV < T s < 200 keV. We have suggested 
that the spectrum of radiation emitted during the breaking 
out of supernova shocks from the stellar envelopes of Blue 
Super Giants and Wolf-Rayet stars, which reach v/c > 0.1 
for reasonable stellar parameters, may include a hard compo- 



nent with photon energies reaching tens or even hundreds of 
keV. This may account for the X -ray outburst associated with 
SN2008D dSoderberg et al.ll2008h . and possibly for other SN- 
associated outbursts with spectra not extending beyond fe w 
100 keV [e.g. XRF060218/SN2006aj dCampanaet alJ2 006)l. 

In this paper we derive exact numerical solutions for the 
steady state structure of RRMS, propagating into a cold up- 
stream plasma of protons and electrons, for shock Lorentz 
factors < 30 and upstream proper densities <C 10 25 cm" 3 . The 
solutions are obtained using a novel iteration method for self- 
consistently solving the energy, momentum and particle con- 
servation equations along with the equation of radiation trans- 
port. We assume that the electrons, positrons and protons may 
be described as a fluid, that the (plasma rest frame) energy 
distribution of positrons and electrons is thermal, and that the 
protons are cold. The v alidity of these assumptions is dis- 
cussed in detail in § 12.31 The Radiation mechanisms that are 
taken into account include Compton scattering, pair produc- 
tion and annihilation and Bremsstrahlung emission and ab- 
sorption. Other radiation mechanisms, e.g. double Compton 
scattering, are shown to have a minor effect on our results. 

The paper is organized as follows. In section §|2]we review 
the physics of RMS, analyze qualitatively the shock structure, 
and motivate the main assumptions. In section §[3] we write 
down the conservation and transport equations that are numer- 
ically solved, in physical and dimensionless forms. In section 
§|4]we present the numerical iteration scheme used to obtain 
the solutions and apply it to several test cases. In section §|5] 
we present the numerical solutions of the shock structure and 
spectrum. In section §[6] we give a simple analytic description 
of the structure of the shock, which reproduces the main re- 
sults of the numerical calculations. In section §|7]we present, 
for completeness, a preliminary detailed numerical solution of 
a non relativistic RMS, and compare it with previously known 
results. In §|8]we summarize the main results and discuss their 
implications. 

Throughout this paper the subscripts u and d are used to de- 
note US and DS values respectively. The term "shock frame" 
refers to the frame at which the shock is at rest (and in which 
the flow is stationary), and the term "rest frame" refers to the 
local rest frame of the plasma, i.e. the frame at which the 
plasma is (locally) at rest, n stands for number density of a 
species of particles, and if not mentioned otherwise refers to 
protons. A summary of the notations repeatedly used in this 
paper appears in appendix § [A] 

2. THE PHYSICS OF RMS 

In this section we discuss the physics of RMS. In § !2.1l we 
define RMS, write down the global requirements that must be 
satisfied by a physical system in order to allow the formation 
of RMS, and derive th e as ymptotic DS conditions. We then 
focus on NR RMS in § 12.21 writing down the assumptions un- 
der which our analysis is carried out, describing the physical 
mechanisms at play, and providing order of magnitude esti- 
mat es fo r the shock width and temperature. Most of the results 
of § l2.2l mav be fou nd in earlier pap ers (Zel'dovich & Raizer 
1966; Weaver 1976; Katz ei1d1l2010h . The physics of RRMS 
is discussed in § 12.31 We highlight the main differences be- 
tween the relativistic and the NR cases, and describe the as- 
sumptions under which the analysis of subsequent sections is 
carried out. 

2. 1 . Introduction to RMS 



2.1.1. Radiation domination 

Consider a steady state shock traveling with velocity cj3 u 
through an infinitely thick, cold plasma of protons and elec- 
trons, with US rest frame density n u . The thermal and radia- 
tion pressures in the asymptotic far DS, which are determined 
by conservation laws and thermal equilibrium, are given by 
2n d T d and obbT^/3 respectively, where n d and T d are the far 
DS proton density and temperature, and orb = tt 2 /15(Hc) 3 is 
the Stefan-Boltzmann energy density coefficient. The radia- 
tion pressure grows much faster than the thermal pressure as 
a function of (3 U , and at high enough /?„ the DS pressure is 
dominated by the radiation. The condition for radiation dom- 
ination is 

obb T d 



» 1n d T d , 
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To obtain Eq. (0, note that at low shock velocities, where the 
radiation pressure is negligible, T d ~ e and n d w 4n„, where 
e fa f3 2 m p c 2 /2 is the kinetic energy per proton in the US. 

2.1.2. Global requirements from a system through which a RMS 

propagates 

In order to sustain a quasi steady state RMS, the system 
which the shock traverses has to be larger than the shock 
width. The width of the deceleration, L^ c , is f3~ l Thomson 
optical depths for NR shocks (see § I2.21 i and, as we show 
in this paper, is ~ 1 Thompson optical depths for relativistic 
shocks. Hence, systems which RMS traverse much satisfy 



L > L dec = (a T nj3) 1 , 



(4) 



Where L is the size of the system, n is the proton density and 
f3c is the shock velocity. The minimum total energy and mass 
of such systems are 

E-^m^nL^-^^SxlO^r'erg (5) 
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respectively. For example, Eq. (0 implies that at ISM typical 
densities, n -C 10 4 , a solar mass rest energy can not drive a 
RMS. In such cases, the shock would be mediated by other 
mechanism e.g. collective plasma processes. 

2.1.3. Far DS conditions 

In the far DS, which is in thermal equilibrium, the condi- 
tions are completely determined by conservation of energy, 
momentum and particle fluxes, 

idFdPd = n u r u /3 u , 

AT 2 d (i d p 1 j = T U I3 U (Y U -T d )n u m p c 2 , 



(4T 2 d f3 2 d + l)p y j = T u f3 u {T u p u -T d f3 d )n u m p c 2 , (7) 

where p ld = l/3asBT d is the far DS radiation pressure, and 
where the plasma pressure in the DS was neglected. Eqs. (J7j 
can be solved for (3 d and T d . In the NR and ultra relativistic 
limits the solution reduces to the expressions 
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respectively, where n u = 10 I5 n„.i5 cm" 3 . Note that the condi- 
tion for a radiation dominated DS, Eq. (|5), can be obtained by 
comparing Tj with e. 

2.2. Non relativistic RMS 

We next focus on NR RMS. By non relativistic shocks we 
refer to shocks in which neither the protons nor the electrons 
move with relativistic bulk or thermal velocities throughout 
the shock. In particular, this implies that the temperature is 
always much smaller than m e c 2 . 

2.2.1. Assumptions 

The discussion below of NR R MS is valid under the follow- 
ing assumptions (Weaverl 1976b . 

• The electron fluid and the ion fluid move together with 
the same velocity. This is justified by the presence of 
collective plasma modes. In the simplest case of pro- 
tons and electrons, an electrostatic field is sufficient to 
couple the fluids. 

• The pressure is dominated by radiation throughout the 
shock transition. This is justified at the end of this sub- 
section. 

• For typical photons, the optical depth is dominated by 
Compton scattering. This, combined with the low ve- 
locity of the flow implies that the diffusion equation can 
be used to approximate the spatial transport of the radi- 
ation. 

• The Compton y parameter is much larger than 1 
throughout the flow, implying that the energy density 
is dominated by a component having a Wien spec- 
trum, and that the electron energy spectrum is close to 
a Maxwellian, with approximately the same tempera- 
ture. The radiation is well described by two parame- 
ters, the temperature and the density of photons n ll& s 
in the Wein-like component. Note, that a large Comp- 
ton y parameter is sufficient to ensure that the electrons 
are strongly coupled to the radiation since the radiation 
dominates the thermal energy density. 

• The main source of photon production is thermal 
bremsstrahlung. 




FIG. 1. — A schematic description of the structure of a fast NR RMS, in 
which the radiation departs from thermal equilibrium. 



2.2.2. The shock transition width 

Physical quantities approach their far DS equilibrium val- 
ues on length scales, which may vary by orders of magni- 
tude for different quantities. In particular, as explained below, 
the transition width of the velocity is determined by Compton 
scattering and occurs on length scales, which may be much 
smaller than the temperature transition width, which is de- 
termined by photon production. A schematic cartoon of the 
vel ocity a nd temperature profiles of NR RMS is shown in 

fig.EZH 

Velocity transition — For NR R MS the width L dec of the ve- 
locity transition region (see fig. 12.2.21 ) is comparable to the 
distance L dt ff ~ (fi u n e tjTY l over which a photon can diffuse 
against the flow before being advected with the flow. To see 
that the velocity transition width can not be larger, note that 
once a proton reaches a point in the shock where the energy 
density is dominated by photons, it experiences an effective 
force 



-r m pC 2 ~ OT$vfi~t ~ cF T n u /3lm p c 2 



implying a deceleration length of 
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The drag estimated in equation ( TTOb is unavoidable due to the 
fact that once the photons dominate the pressure, they cannot 
drift with the protons, as this will imply a radiation energy 
flux greater than the total energy flux. 

Thermalization length — The region of the shock profile over 
which the temperature changes before it reaches Td can be 
extended to distances that are much larger than L<#//. To see 
this, consider the length scale that is required to generate the 
density of photons of energy ~ T d in the DS, determined by 
thermal equilibrium, n 7 , eq ps p 7t d/Tj, 



f3c 



QfjeS 



(12) 



where Q-y >e s is the effective generation rate of photons of en- 
ergy 3T d . We use here the term "effective generation rate" 
due to the following important point. Photons that are pro- 
duced at energies <C T d may still be counted as contributing 
to the production of photons at T d , since they may be upscat- 
tered by inverse-Compton collisions with the hot electrons to 



energy ~ Td on a time scale shorter than that of the passage 
of the flow through the thermalization length, Lj / (3dC. The 
Bremsstrahlung effective photon generation rate is given by 



Q-yfiS = a e n p n e a T c 



-Aeffgeff) 



(13) 



where g e ff is the Gaunt factor, A e ff ~ \og[T / (hv m i n )] and v m i n 
is the lowest frequency of photons emitted by the plasma 
which may be upscattered to 37^ prior to being absorbed (ab- 
sorption is dominated by Bremsstrahlung self absorption for 
far DS values). 
The resulting thermalization length is 



L T (3 u n u a T 



1 



(14) 



(15) 



lOOaeAeffgefE y/m e c 2 T d m p c 2 
This implies that for high shock velocities, 

U > 0.07«;/ 30 (A eff £ eff ) 4 / 15 , 

the length required to produce the downstream photon density 
is much larger than the deceleration scale. For lower shock 
velocities, thermal equilibrium is approximately maintained 
throughout the shock. 

2.2.3. Description of the shock structure 

An analytic expression for the velocity, density and pres- 
sure profile can be found under the diffusion approximation 
(e.g. IWeaverlll976h . In particular, the velocity (3c at a give 
position x along the shock satisfies 
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(16) 



The shape of the temperature profile is largely determined 
by the photon production in one diffusion length into the DS 
(the first (3j l optical depths of the downstream region, hence- 
forth the immediate DS). In this region the photons mediating 
the shock are produced. If a photon density of ~ obbT] /3 
is produced {Lj < Ldiff), the flow will stay close to thermal 
equilibrium, and the temperature profile, which can be ex- 
tracted directly from the analytic pressure profile, essentially 
follows the velocity profile. Otherwise, when Lj 3> Ldiff, the 
velocity transition of the shock ends without reaching the far 
DS equilibrium temperature. The radiation pressure n lfi gT 
reaches its DS value as soon as the velocity is close to the 
DS velocity. Down stream of this region, the density of pho- 
tons « 7i eff grows with distance as more and more photons are 
being generated and advected with the flow and saturates at 
the equilibrium black body photon density sa obbT] /3 . Ac- 
cordingly, T is decreasing throughout the downstream. In this 
case, we can broadly divide the shock structure into four sep- 
arate regions. 

1. Near upstream: A few diffusion lengths, (f3 s a T n u )~\ 
upstream of the deceleration region. In this region, 
characterized by velocities that are close to the up- 
stream velocity, j3 « /3 U , and temperatures T ^> T u , the 
temperature changes from T u to ~ T s . It ends when the 
fractional velocity decrease becomes significant. 

2. Deceleration region: A (f3 s o-Tti u )~ l wide region where 
the velocity changes from j3 u to fid and the temperature 
is roughly constant, T ~T S . 



3. Immediate downstream: Roughly a diffusion length, 
(j3 s o-T n u)~ l , downstream of the deceleration region. In 
this region, characterized by velocities close to the 
downstream velocity, f3 sa f3 d , and temperature T ~ T s , 
the photons that stop the incoming plasma are gener- 
ated. Upstream of this region j3 > j3 c { and the photon 
generation rate is negligible. Photons that are generated 
downstream of this region are not able to propagate up 
to the transition region. To estimate the temperature 
value in the immediate DS, T s , the number of photons 
produced in the immediate DS by Bremsstrahlung and 
up-scattered by inverse Compton should be equated to 
the number of photons required to carry the pressure 
at that point. The production rate, given by Eq. (fT3l . 
combined with diffusion and conservation laws, leads 
to the following estimate of the immediate DS temper- 
ature for NR RMS ( Kat z et al.H20Tob 
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where A e ff,i = 10A e ff. This resu l t is in agreement with 
the numerical results of lWeaverj (119761) . 

4. Intermediate downstream: The region in the down- 
stream where most of the far downstream photons are 
generated and T changes from T s to T c {. This region 
has a width Lj given by Eq. ( TBI , much grater than 
(/SjOrW,,) -1 . Thus, diffusion within this region can be 
neglected. The temperature profile is expected to fol- 
low T oc x~ 2 . To see this, note that the photon den- 
sity at a distance x from the shock is proportional to 
the integral of the photon generation, « 7l eg oc T~ l l 2 x. 
Since the photon pressure equals the downstream pres- 
sure, we have « 7j eff oc T~ l and T oc x~ 2 (this is valid 
for a constant value of A e ffg e ff an d is somewhat shal- 
lower in reality). Using this dependence of the temper- 
ature on distance, the thermalization length can shown 
to be related to the deceleration length by Ljfi s n u aj ~ 



-1/2 
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2.2.4. Scaling of the profile with density 

The velocity, density and pressure profiles of RMS (as a 
function of optical depth r = ajxriu) are independent of the 
upstream density. The scaling of the temperature depends 
on whether or not thermal equilibrium is sustained. In case 
it is (if Lj <C Ldiff), the temperature scales with density as 

T oc nj . Alternatively, when Lj 3> Ldiff, in the shock re- 
gions where the temperature is much higher than its equilib- 
rium value, the temperature profile does not scale with den- 
sity, T oc n° u , see for example eq. ( fTTI i. To see this, note that in 
these regions Bremsstrahlung absorption is negligible, while 
Compton scattering and Bremsstrahlung emission are both 
two body processes that scale similarly with n u . The conser- 
vation and radiation transfer equations are invariant under the 
scaling of the radiation intensity, densities and length scales 
across the shock by n\ ,n\ and n~ l r espe ctively. This scaling 
is shown explicitly later, in section § 13.11 Bremsstrahlung ab- 
sorption may still be important at low frequencies and affect 



the structure through logarithmic corrections to the effective 
photon production rate. In the far downstream, where T ap- 
proaches its equilibrium value, absorption will no longer be 
negligible. 

2.2.5. From radiation domination to radiation mediation 

For shocks satisfying Eq. (0, the far DS pressure is dom- 
inated by radiation. It is not a priory trivial that the pressure 
is dominated by radiation in the velocity transition region of 
such shocks since photons that are generated in the DS are 
able to diffuse upstream over a finite distance only. We next 
illustrate that under a wide range of conditions, the radiation 
does indeed dominate the pressure in the velocity transition 
region. 

Consider a hypothetical shock, having a DS energy den- 
sity dominated by radiation, in which the velocity transition 
is mediated by some mechanism other than radiation. In the 
absence of radiation, the temperature immediately behind the 
velocity transition would be T ~ e = 0.5(3 2 m p c 2 . Photons gen- 
erated in this region can diffuse upstream to a characteris- 
tic distance of Ldiff > {(3n e <jT)~ l ■ Under these assumptions, 
the energy in photons that are produced by Bremsstrahlung is 
much larger than the available thermal energy, 
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where Qs r is the photon production rate by Bremsstrahlung 
at energy ~ T, and where we used /3 = /3„/4, appropriate for 
NR shocks which are not radiation mediated and e ~ T ^C 
m e c 2 . Note, that by definition, the condition (TT8T > is roughly 
equivalent to demanding that the temperature T s in ( fT7T i be 
smaller than e. 

This implies that a shock with a radiation dominated DS and 
negligible radiation in the velocity transition region cannot 
exist if the transition region is smaller than Luff. Once the 
pressure is dominated by photons, they will also mediate the 
shock. 

2.3. Relativistic RMS 
2.3.1. Assumptions 

Throughout this paper we make the following assumptions 
forPvRMS: 

1 . The electrons, positrons and ions move as a single fluid 
with the same velocity. This is motivated below by the 
presence of collective plasma instabilities. 

2. The electron and positron velocity distributions in the 
rest frame are approximately thermal. This assumption 
is justified by the intense radiation field interacting with 
the electrons and positrons, which quickly eliminates 
large deviations from the mean velocity. 

3. The ions have a negligible contribution to the pressure. 

4. The radiation mechanisms dominating the shock are 
Compton scattering, bremsstrahlung emission and ab- 
sorption and two photon pair production and pair anni- 
hilation. 

The assumption of a single plasma velocity is motivated by 
the fact that the plasma time (t p i) is much shorter than the 



mean time between Compton scatterings (t c ) of an electron, 
allowing for collective plasma processes to isotropize the ve- 
locities of the particles. Indeed, the ratio of these timescales, 
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where n e = I0 l9 n e 19 cm" 3 , is much lower than unity given that 
tij/tie is not very large, see § 12.3.31 

The second assumption we make regarding the plasma - the 
existence of an effective temperature, is somewhat more sub- 
tle. In principle, the electrons and positrons can have a general 
distribution function. However, most of the shock is char- 
acterized by a strong dominance of radiation energy density 
over particle thermal energy density. This leads to the elec- 
trons and positrons being "held" in momentum space by the 
radiation, since each scattering changes the energy of the elec- 
tron considerably, if it departs significantly from the average 
photon energy. The only way to maintain a very non-thermal 
electron spectrum is by having a radiation spectrum which is 
not dominated by a typical photon energy, e.g. a power law. 
Our numerical results show that the radiation energy density 
is dominated by photons of limited energy range in the rest 
frame of the plasma, and that when a high energy photon 
tail appears, the photons populating this tail have a very low 
cross section for interaction with electrons or other photons. 
This supports the assumption that the energy distribution of 
electrons and positrons may be characterized by some typical 
"thermal" energy, greatly simplifying the calculations. 

A note is in place here regarding Coulomb collisions. The 
effective cross section for Coulomb collisions of electrons on 
protons is a ~ e 4 /ej, where £k is the electron kinetic energy. 
When the energy of the electron is of the order of m e c 2 , the 
cross section is similar to oy, the Thomson cross section. This 
implies that Coulomb collisions in RRMS play a marginal role 
in equilibrating the motion of particles in the plasma, as the 
photon density inside the shock is typically of the order of the 
electron density. Unlike plasma instabilities, this process can 
not account for the equilibration of the distribution function 
of the particles. At low energies, i.e. NR RMS, Coulomb 
collisions may become dominant (see Weaver 1976) due to a 
much larger effective cross section. 

2.3.2. Velocity and temperature transition regions' widths 

Velocity transition — The line of arguments presented in § 12.21 
for estimating the velocity transition width can not be di- 
rectly extended to relativistic shocks since KN corrections to 
the Compton scattering cross section depend on the a priori 
unknown photon frequency and plasma temperature, which 
vary throughout the transition region. Note, that as expa- 
lined above, pair production and relativistic corrections to the 
cross sections become important already at non relativistic up- 
stream energies e = f3 2 m p c 2 /2 ~ 100 MeV, since the temper- 
ature of the plasma within the deceleration region reaches a 
considerable fraction of m e c 2 for this value of e. The pro- 
duction of pairs also changes the simple estimate, since it 
changes both the scatterers' number density and the shock op- 
tical depth. Finally, an additional complication is introduced 
by the strong dependence of the scattering mean free path on 
the photon's direction of propagation, expected due to the rel- 
ativistic velocity of the plasma. 



Thermalization length — The thermalization length can be es- 
timated in a way similar to the NR case, since the width of 
the temperature transition is much larger than the deceleration 
width [see Eq. (1151) 1. and the scale is set by the lowest tem- 
perature, i.e. r</, which is non relativistic. The thermalization 
then takes place over 



1 



100a c A e ffgeff ^Jm e c 2 T d m p c 2 



(20) 



Thomson optical depths. Since s = (T„- l)m p c 2 and Tj is 
many orders of magnitude smaller, this width is always very 
large in terms of Thomson optical depths. 

2.3.3. Immediate DS 

We next give a rough estimate of the average temperature 
in the first f ew optical depths of the immediate DS of RRMS 
dKatzet al] 2010). The assumption we use is that the electron- 
positron pairs and the radiation are in Compton Pair Equi- 
librium (CPE). This assumption is valid since the velocity is 
< c/3, and since the y parameter arising from mildly relativis- 
tic temperatures is large, as shown below. The numerical cal- 
culations are not based on t his assumption, and its self consis- 
tency is discussed in § 16.11 Following the NR RMS analysis, 
T s is estimated by equating the number density of photons pro- 
duced by Bremsstrahlung and by inverse Compton emission 
of thermal pairs with the number density of photons needed 
to carry the energy flux at the end of the deceleration region. 

Assuming that the number density of pairs is much larger 
than that of protons, and neglecting Double Compton emis- 
sion, the ratio of photon to electron-positron number densities 
may be written as 



= -a e A eS g e f fM (T)P d , 

rii 5 



(21) 



where » 7>e g is the density of photons in the Wein-like compo- 
nent (see § 12.2. U and the free-free emission is written in the 
form 

Q 7fi s= a e a T cnfA eB g ef f 7ie i(f). (22) 

Here geff.rei is the total Gaunt factor [defined by Eq. ( f22l l 
including all lepton-lepton Bremsstrahlung emission. For 
10 < A e ff < 20 and 60 keV < T < m e c 2 , the approximation 

feff,rei ~ A e ff/2 (23) 

agrees with the results of SvenssonJ (119841) to an accuracy of 
better than 25%. At these high temperatures, the Compton y 
parameter is large and radiative Compton emission is negligi- 
ble. Substituting Eq. ( 1231 in Eq. (fJTJ we find 

2 

(24) 
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«7£» _ 9 r / iV eff \ „ -2 



In the regime 200 keV < T < m e c 2 , pair production equilib- 
rium is approximately given by 



rij/ni « Q.5m e c /T. 



(25) 



Comparing equations d24b and ( |25] >. we see that if 
T > 200 keV there would be too many photons generated per 
lepton. Much lower temperatures lead to insufficient photon 
production, as can be deduced from the NR case. We con- 
clude that for relativistic shocks, 



r 5 -200keV. 



(26) 



The weak dependence of the immediate DS temperature on 
parameters is due to the rapid increase of pair density with T 
at T ~ m e c 2 . 

Subsonic region — The US flow is "super-sonic", in the sense 
that the plasma velocity, f3 u c, is larger than the plasma speed 
of sound, (3 ss c. The production of a large number of pairs in 
the immediate DS, n+/n p S> 1, and the heating of the plasma 
at this region to relativistic temperatures, 7^ ~ QAm e c 2 , im- 
plies a "sub-sonic" flow, /3 SS > (3, in the immediate DS. The 
large number of pairs implies that the average plasma particle 
mass is close to m e , for which the temperature, T s ~ 0Am e c 2 , 
gives a speed of sound which is close to it's highly relativistic 
value, j3 ss = l/v3 (see appendix §|C]for a detailed calculation 
of the speed of sound). j3 ss = I/a/3 is larger than the plasma 
velocity in the immediate DS, which is close to its far DS 
value, p d < 1/3. 

Note, that we are referring here to the plasma speed of 
sound neglecting the (dominant) contribution of the radiation 
to the pressure. This speed of sound describes the propaga- 
tion of (small) disturbances in the plasma on length (time) 
scales which are short compared to the mean free path (ti me) 
for electron-photon collisions. As explained in § 12.3.11 see 
eq. ( fT9l , collective plasma modes are expected to lead to a 
fluid like behavior of the plasma on length and time scales 
much shorter than the electron-photon collision mean free 
path. 

In the far DS, the flow becomes super-sonic again, f3 ss < /?</. 
This implies that for relativistic shocks the flow crosses two 
sonic points, accompanied by singularities of the differential 
conservation equations [Eqs. ( 1281 ). d29|)l: 

• At the first sonic point, the flow changes from super- 
sonic to subsonic. This is a hydrodynamically unsta- 
ble point which results in a hydrodynamic shock. A 
steady state hydrodynamic flow can not smoothly cross 
a sonic point going from supersonic to subsonic veloc- 
ities because downstream of the sonic point, upstream 
going characteristics conv erge to the sonic point (e.g. 
Zel'dovich & Raizer 1966), infinitely steepening a con- 
tinuous profile at that point and resulting in a shock. 
We show below that while most of the deceleration of 
the plasma is continuous, a (sub-)shock across which 
the velocity jump is small, S(Tf3) ^0.1 (see fig.|7]), is 
indeed required to exist at the end of the deceleration 
region. This sub-shock must be mediated by the same 
processes that are assumed to isotropize the particles' 
velocities in the fluid rest frame on a scale much shorter 
than the radiation mean free path [e.g. plasma instabil- 
ities, see eq. $1%]. 

• At the second sonic point the flow passes from a sub- 
sonic to a supersonic region. This is a stable point 
which has no special significance, and is simply part 
of the thermalization tail of the shock. 

2.3.4. Structure 

The structure of RRMS differs from that of NR RMS. The 
main differences are: 

• The deceleration length is much larger than the naive 
estimate: The length, measured in Thomson optical 
depths of e~ e + , grows with the upstream Lorentz fac- 
tor r„ in a manner faster than linear (for NR RMS it is 

-AT 1 ); 



• Pair production has a significant contribution to the de- 
celeration of the plasma. 

• As explained above, a "hydrodynamic" sub-shock (pos- 
sibly mediated by plasma instabilities) across which the 
velocity jump is small, 5(T/3) ~ 0.1, is required to exist 
at the end of the deceleration region; 

• The radiation is highly anisotropic, and exhibits a high 
energy tail with a typical cutoff energy of ~ T 2 m e c 2 . 

2.3.5. From Radiation domination to radiation mediation 



Expanding the reasoning given in § |2.2.5| for NR RMS, we 
argue here that relativistic shocks which are radiation domi- 
nated, i.e. in which the DS energy density is dominated by 
radiation, must also be radiation mediated. Let us assume 
the contrary, i.e. that the energy density in the deceleration 
region is not dominated by radiation and that deceleration is 
therefore not mediated by radiation. In this case, the plasma 
reaches a temperature ~ T u m p c 2 at the end of the deceler- 
ation region, and then gradually thermalizes as it flows fur- 
ther into the DS. Since Tj < m e c 2 (see eq.|9|, let us consider 
the point in the downstream where the temperature reaches 
T ~ 10 MeV. Since the velocity at this point already reached 
its DS value, ft ~ 0j < 1/3, photon transport is well described 
in this region by the diffusion approximation, with diffusion 
length Ldi/f ~ 3(n e aT)~ 1 . An electron crossing this diffusion 
length produces a large number of ~ 10 MeV photons, 

OC e gff,relKff 



»7,10 



n e f3c 



Q 



7>e// ' 
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10 8ffM(10MeV)A eff 
100 



(27) 



(using conservative estimates for the Gaunt factor and the log- 
arithmic correction). This ratio is much larger than its CPE 
value n 7 .eff ~ n e , expected at T ^> m e c 2 . For such a high ra- 
tio of photons to electrons, a photon will produce a pair on 
another photon on a time scale much shorter than its scatter- 
ing time scale. Such a deviation from equilibrium on a length 
scale ~ 3 scattering optical depths is not self consistent. We 
conclude therefore that the temperature can not significantly 
exceed m e c 2 at the point where the deceleration is complete. 
This implies, in turn, that most of the energy at the end of the 
deceleration must be carried out by radiation. 

3. RMS EQUATIONS AND BOUNDARY CONDITIONS 

bf structural change in this section. Intro added. In this 
section we write down the equations of RMS that a re num er- 
ic ally solved based on assumptions 1-4 given in § 12.3.11 In 
§ 13.11 we write down the hydrodynamic and radiation trans- 
fer equations in physical and dimens ionle ss form and define 
the variables that are solved for. In § !3.2l we provide expres- 
sions for radiation scattering (Compton), production and ab- 
sorption (Bremsstrahlung) and pair production and annihila- 
tion. A sum mary of all the equations in dimensionless form 
is giv en in § 13.2.41 The boundary conditions are described in 

3.1. Hydrodynamic and radiation transfer equations 

The equations governing the structure of a steady planar 
shock propagating along the z direction are 



— T 0z -Q 



(28) 



dz s h 

r„fl, 

n P - n PM po J 

d(Tpn + ) _ Q + 
dz s h c 



(29) 
(30) 
(31) 



H>sh' 



dluJjUshi 



dz. 



■sh 



' llshiHslnVsi^-IvJPshlXshiHslnVsh). (32) 



The first 3 eqs. describe the conservation of energy, momen- 
tum and proton number. The forth eq. describes the pro- 
duction and annihilation of positrons, and the fifth eq. de- 
scribe the transport of photons. z s h is the shock frame distance 
along the shock propagation direction, j3c is the plasma ve- 
locity in the shock frame, V = 1/ yl — ft 2 is the corresponding 
Lorentz factor, n p is the proper proton density, n+ is the proper 
positron density, I Vth (n s h) is the shock frame specific inten- 
sity at (shock frame) frequency v s i, and direction /i,/, = cos 6 s f, 

(9 s i, is the azimuthal angle with respect to z), and T^f is the 
shock frame energy-momentum tensor. Q + = dn + /dt is the 
net positron production rate (production minus annihilation), 
which is frame independent. 77 and \ stre the emissivity and 
absorption coefficients, respectively, and are functions of the 
plas ma parameters and of the local radiation field described 
in § 13.21 We use {v,I v ,fi,rj,x} to denote quantities measured 
in the plasma rest frame, and add a subscript "sh" to denote 
values of these quantities measured in the shock frame. 

The energy and momentum are carried by the plasma and 
the radiation, 

1 sh ~ l sh,pl T l sh,radi ^ JJ> 

where the subscripts pi and rad refer to the plasma and ra- 
diation contributions respectively. The radiation part of T s "f , 

Pradjh, is given by 



sh.rad 



■■ F rad , sh and 1%. : 



F, 



ad ' 



rad ' 



dQfidvIvip), 



dQfi dvl v (\i) 



(34) 



(35) 



(see appendix D for rules of transformation between rest 
frame and shock frame measured quantities). As mentioned 
in the introduction, we assume that the protons, electrons 
and positrons may be described as a fluid of single veloc- 
ity cp(zsh), that the energy distribution of the electrons and 
positrons is thermal, with temperature T(z s h), and that the pro- 
tons are cold. Under these assumptions, 



and 



T° P l sh = T 2 P{e pl+ P pl ) 



T'f sh =P P i + r 2 p 2 (e pl+ P l 



Pi) 



(36) 



(37) 



where the proper energy density e p i and pressure P p \ are given 

by 



e p i = n p m p c +(n e + n+)m e c +-f(T)(n e + n + )T, (38) 



and 



P p i = (n e + n+)T 



(39) 



where 




FIG. 2. — A comparison of the exact value of f(T) (solid line), calculated 
num eric ally for a Maxwellian distribution, and the approximation given by 
Eq. go) (dashed line). 

(note, that we neglected the thermal pressure of the cold pro- 
tons). f(T) is dimensionless and is approximated by the fol- 
lowing interpolation between the NR (/ = 1) and relativistic 
if = 2) values, 



/(r) = itanhf ln(r ^ c2) + - 3 V^. 
J ' 2 \ 1.93 / 2 



(40) 



This approximation describes the equation of state of 
Maxwell-Boltzmann distributed plasmas to an accuracy bet- 
ter than ~2x 10~ 3 for all temperatures, as shown in fig. [2] 



3.1.1. Dimensionless equations 
We define the following dimensionless quantities: 
T 





m c c z 
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hv 
m e c L 




x+ 


= n+/n p , 




Zsh 


= r u n u aTZ S h, 




dr* 


= T(1 + I3)(n e 


+ n+)o- T dz s h 


i 


I 
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With these definitions, and using the explicit forms of T , 
derived above, the energy and momentum conservation equa 
tions take the form 



1 « L m,p 

+ 2nF rad , sh = 1 + — , 



YB f m e 

=-g- l + (l+2x + )^ 
Y u Pu I m p 

1 - m e 

+ —%KP ra d,sh = 1 + ■ 



\+T[l + -f(T) 
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(44) 



(45) 



are the scaled energy and momentum fluxes of the radiation 
field. 
The transfer equation, eq. (l32l , takes the form 

Ush—^ — — =Vsh(Vsh 1 £ , sh)-h< h (p> s h)Xsh{p,sh,v sh ). (46) 

dT* 

The emissivity and absorption coefficients are the sum of the 
contributions due to the various processes considered 



fltotiV, V) = Z_. VprociH, *>), 



Xtot 



(M^) = ^2xprocQi,0)- 



(47) 



(48) 



The transformation relations for the scaled emissivity and ab- 
sorption are 



rj = 



>1 



Y( 1 + /3)a T (n e + n+) m p Y 2 u /3 u n u hc ' 
X 



x- 



Y(\+P)o- T {n e + n + Y 



(49) 



(50) 



Finally, the equation describing the evolution of pair density 
may be written as 



dx+ 

-t ± = Q+, 

dr* 
where the scaled rate of pair production is 

Q + 



Q + 



(51) 



(52) 



Y 2 B( 1 + f3)n p (n e + n+)a T c ' 
We describe next the various radiative processes included. 

3.2. Radiation mechanisms 

The radiative processes we take into account are Compton 
scattering, Bremsstrahlung emission and absorption and two 
photon pair production and annihilation. Other processes, 
which we neglect, do not modify the results significantly. 
The leading corrections are due to double Compton scattering 
(7 + e — > 2j + e), three photon pair annihilation (e + e~ —> 3j) 
and pair production on nuclei. Other processes, such as muon 
and pion pair production and synchrotron emission, are less 
significant. 

3.2.1. Compton scattering 
The contribution of Compton scattering to rj and \ is 

/j 
dn ' dv 'd~7akv ( V ' ,Q ' ~* "'") Il " (n ' ): 

(53) 

Xs^v) = {n e +n+)o- c {v,T)y. (54) 



where the total cross section 

/ 

dv'dVL 



o c (y,T) = I dtl'di/ -^- (v,n-n/,Sl'), (55) 



may be written as 

cT c (OJ)=Jdttdv'-^^ (i/,!)4 v',n') 
i + Cf 2C(i + C) 



-o~t~ 



C 3 I 1+2C 



ln(l+20 1+3C 



-ln(l+2C) 



2C (1 + 2Q 2 



(56) 



Here, C = 0(1+2T) [see e.g. iRvbicki & Lightmanl d 19791) 1. 
The normalized emissivity and absorption are 



[T(l-j3/j, sh )] fj s ,sh(O s h^sh) = Vs(0, 51) = 



da 



1 



r(i+/3) 



x dn'do'-—j-{o',n'^ o^h^'), (57) 



^wft 



[r(l - ^a)]- 1 XsAOsh) = Xs(0) = r(1 1 + a c (0, f ), (58) 

where er = cr/aY and the transformations between the shock 
frame and plasma rest frame values of v, ti are given in ap- 
pendix § |D] 

Since using the exact form of the differential cross sec- 
tion for Compton scattering greatly increases the compu- 
tational resources demands, we use instead an approxima- 
tion described in appendix § [B] In particular, we assume 
isotropic scattering in the rest frame of the plasma, i.e. 
da s (v, 51 — > v' , 51') independent of 51'. 

3.2.2. Pair production and annihilation 

Pair annihilation — The photon emission arising from anni- 
hilation of pairs has the form 

1 . hvn e n + a T cf±(v,T)r±(T) 
T\ v =—n v hv= , (59) 

4-TT 4-TT 

where r± is a dimensionless function of T accounting for the 
rate of annihilation and f± is the spectral distribution of the 
photons, where 



f±(vj)dv=l. 



(60) 



The approximati on we use for f± is based on the analysis of 
Zdziarski (1980), who fits an analytic function to the results 
of Mont e Carlo calculatio ns. For the annihilation rate we use, 
based on Svensson ( 1982b, 



r ± (T)=- 



1 + 



27/ 2 



ln(2ry £ r+1.3) 



(61) 



where tje = <? 7£ ~ 0.5616, and je ~ 0.5772 is Euler's con- 
stant. The normalized emissivity is given, based on Eq. 
by 



fl±(0, 51) = [T(l - /3/Ish)] fl±,sh(Vshi^lsh) 

_ (x+ + l)x+0f±(0,f)r ± (T) m e 
~ 47r(2x + +l)r„r 2 /3(l + /3)m/ 

The annihilations rate in Eq. ( l3TT i is simply 



Q = --n e n + <j T cr±(T), 



(62) 



(63) 



and the scaled contribution to Eq. ( BTT l is 

x + (x + +l)r±(T) 



Q + 



2r 2 /3(l+/3)(l+2x + )' 



(64) 



Pair production — The two photon pair production contribu- 
tion to the absorption in the transfer equation is 



Xu, T y(l^)= / cr 77 (^,^',M ; 51')x 



I ^——(l-cose l )<d[vv'(l-cos6x)-2vl}dQ,'dv\ (65) 
chv' ' 

where 0\ is the angle between /i and 51'. The scaled absorption 
can be written as 

X ^ (P)= (1 + P)(2X ++ 1) J a ^ V '"'" )X 

lv \ l-cos9i)Q[00'(l-cos6i)-2]dn'd0' (66) 
v' 

(X should be calculated at the same frame f or which I' is 
given) . For the cross section we use [e.g. Padmanabhan 
(2000)] 



cr 77 (s) = 



3 a T 
87 



where 



2+ ^_ _L Wy,_ 1+ I )(,-! 



s= — hvhv'(l — nix') 



1/2' 



,(67) 



(68) 



is the center of momentum energy squared. To shorten 
the computing time we integrate over <j>' assuming that <r 77 
changes slowly with <p' and that Q[vv'(l-co?,6\)-2v 2 p \ has 
the same value for most </>' values, obtaining approximately 

< l-COS(6>i) >0= 1-flLl'. 

To find the positron production rate Q+ we use the rate of 
photon loss to this process, 

2+ = ~2«7 = 2 / \ — Xvwil^dvdQ. (69) 

The scaling of the production rate follows, 



T u m p f Io(ijl) 

Q+ = -= / —^Xtaid^dOdVL. 

2m e J v 

3.2.3. Bremsstrahlung 



(70) 



Bremsstrahlung emission includes contributions from e~p 
and e + p encounters, as well as from e~e~,e + e + and e~e + en- 
counters, which become important sources of photon produc- 
tion at high temp eratures. The emission can be expressed by 



tion at nign tempe 
dSvenssorJ 19821) 



1 FT e~ hv l T 

tiryfAQjO) = —z\ — a e a T m e ' c n t — -^— A 

7T Z V 7T y/Tv 



//) 



(71) 



where a e is the fine structure constant, and 
\ ff (x + ,T) = (l+x + )\ ep + [x 2 + + (l+x + ) 2 ] X ee +x + (l+x + )X + 



(72) 
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is a numerical factor accounting for the presence of electron- 
positron pairs and for relativistic corrections at high tempera- 
tur e. We use a prescription for bremsstrahlung emission based 
on Skibo et al. ( 1995) (note that there is an errata correction 
to this paper), which gives a general fit for the Gaunt factor 
as a function of temperature, positron density and the emit- 
ted frequency. The transformation between the different nota- 
tions is X/f = -rj^gs, where g s is the Gaunt factor as given in 

lSkiboetalJ(fl99l . 

The emissivity resulting from Eq. ( 1711 is 



r lffA^) = hvn 1 jf{Vl,v)-- 



1/2 2 2 e 

a^ajnij c n i — 



-hv/T 



V/- 



7T V 7T Vr 

The normalized emissivity then reads 

[T(l - (lUsh)] 2 VffM^h, Qsh) = *?//(£, ft) ■■ 

a e m e /m p j2 e~ u ' T 



(73) 



(74) 



Minimal v — Coulomb screening suppresses bremsstrahlung 
emission at impact parameters larger than the Debye length 
Ad = y/T/4-ne 2 (n e + n+), implying a low energy cutoff for 
bremsstrahlung emission (Weaver, 1976b) 



7eJhPe,th 



A, 



he . 



(75) 



D 



where j et f, is the Lorentz factor associated with the random 
("thermal") motion of the electrons, and /3 e ,th is the associated 
velocity (in units of c). Setting 7^,/, ~ 1 + 3T/m e c 2 we get for 
the non relativistic case (T <C m e c 2 ) 



2.87 x 1(T 6 n'(g(l +2x + ) 1/2 KeV, 



(76) 



and for the relativistic case (T 3> m e c 2 ) 



€sc,rel 



( T A 3,/2 
9.12 x 10- 10 «, I f 5 (l+2x + ) 1 / 2 f —J KeV, (77) 



where n,- = nt t \$\ 

We note that since our calculation explicitly describes up- 
scattering and bremsstrahlung self absorption, there is no need 
to introduce (as was done, for example, in Weaver 1976) a 
cutoff to the Bremsstrahlung emission at low frequencies, for 
which the flow dynamical time scale or the self absorption 
time scale are shorter than the time required for a low energy 
photon to be upscattered to T. 

Bremsstrahlung self absorption. Using Kirchhoff 's law and 

the calculated value of r\ v jf in the rest frame of the plasma 

we have 

v _ VvJf 
Xvjf ■ 



where 



BAT)- 



BAT) 



2hv i 



[cm" 1 ], 



(78) 
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is Plank's spectrum. 
The normalized Plank spectrum is 



B f ,= 



1 



h?m p Tlp u n ue >lT-\ 



(79) 
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and the normalized absorption is 

[IX 1 - /3^/Or 1 XffM^h) = Xff(v) ~ "' ' ' ' 

a e h 3 \ (ff) 



BAT) 

n, (l-e- p / f 



V2ir 5 / 2 m 3 e c 3 T(l +/3)(1 +2x + ) 



(81) 



3.2.4. Summary 

To summarize: we use equations d42b . d43l , d46b and 
( BTT l. to determine the variables T(r*), /3(r*), x+(t*) and 
UdXl^shXT*)- The contributions of the radiative processes to 
the transfer equation [eq. j46l )l are given by eqs. d5"7| ), d5"8"1 ), 
J62l , j66l ), j74] > and (f8Tb . The contributions to Q in the 
positron fraction equation [eq. ( 15X1 )1 are given by eqs. j64l i 
and <|70]). 

3.3. Boundary conditions 

We obtained solutions of the equations given above over 
a finite optical depth range around the shock transition, that 
satisfies the following requirements: 

• The solution includes a subsonic region downstream 
of a supersonic region with continuous radiation field 
Iv sll {pL s h) and positron flux across the sub-shock separat- 
ing the two regions (see § 12.3.31 ): 

• The radiation momentum flux in the last several photon 
mean free paths away from the shock transition in the 
US region is negligible compared to the far US electron 
momentum flux; 

• The width of the subsonic region is sufficiently large 
compared to the photon mean free path, such that the 
solution is insensitive to the precise boundary condi- 
tions that are applied at the DS edge, while remaining 
short enough as to avoid reaching the second supersonic 
region which exists DS of the subsonic region. 

The boundary conditions in the far upstream are I^Xl^sh > 
0, Zsh = -00) = 0, i.e. no incident radiation at the upstream 
(In practice we use an effective "reflector" in the US end of 
the calculation, to avoid numerical fluctuations and shorten 
the iteration time. It does not affect the shock structure). In 
addition, the positron number is taken as at the US boundary. 

The boundary conditions at the far downstream are given by 
thermal equilibrium. Since the calculation does not reach the 
far DS, we use a boundary condition in the DS which corre- 
sponds to isotropy of the radiation field in the rest frame of the 
far DS. This is done by equating the intensity and spectrum of 
US going radiation at the DS boundary to that of the DS going 
radiation. For numerical reasons, we multiply the reflected ra- 
diation by a factor which is close to unity, this has a negligible 
effect on the shock structure. In addition, we impose an upper 
limit on the photon energy of the reflected radiation, typically 
3m e c 2 . The physical reasoning for this upper limit is that high 
energy photons that cross this point in the DS either scatter 
and lose most of their energy (as T <C 1 at that point and fur- 
ther away), or more likely, produce an e + e~ pair that is swept 
DS with the flow. 

4. THE NUMERICAL METHOD 

We briefly describe below the numerical method we use for 
solving the equations. 
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4. 1 . Iteration scheme 
We start with an initial guess for the shock profile, S° = 
{f o (T*),0 o (T*),xl(T*),I^ sh (n S h,T*)}, and modify it iteratively 
until a solution of the equations is obtained. The iterations are 
performed as follows: 

1. Compute fj(S") and x(5") using the profile S"; 

2. Integrate directly the transfer eq., eq. d46i l. using i](S") 
and x(S"), to obtain T^O^r.); 

3. Use eqs. d42b . d43l and ( BTT i with the new radiation field, 
^ +1 (Ms/i; t *)' t0 obtain the new profile 5" +1 ; 

Usage of "partial iterations", where S" +l is replaced with a 
weighted average of S" and S" +l , was required in order to 
achieve convergence and stability. 

At any given r*, the energy and momentum conservation 
equations, Eqs. ( l42l . ( |43l , have a supersonic and a subsonic 
solution for (3 an d T g iven x+, I& A (fish)- The position of the 
sub-shock (see § I2.3.31 l was set to r* = 0, upstream of which 
the supersonic solution was chosen and downstream of which, 
the subsonic solution was chosen. 

We significantly reduced the computational time of the cal- 
culation, by separating the spatial grid into two regions, and 
preforming the above iterations on each. The downstream go- 
ing photons on the downstream boundary of the first region 
were used as a boundary condition for the second region and 
vice versa. We preformed macro iterations in which we up- 
dated these boundary conditions until a self consistent profile 
was obtained across the border between the regions. 

4.2. Discretization 
We use a discrete approximation of Ish,v, h ((i>sh)> 

hh,v sh i/J-sh) = / J Ishjjffl {fsh , V s hs , V s h,i+\)fr\{P'sh , l^sh.j , fl>sh,j+l ) , 

^ (82) 

where fn(.x,Xi,X2) = ®(x—X\)Q(x2—x) is the top hat func- 
tion and 9 is the step function. The distribution of v s u : i 
is logarithmic in the range v min to v max . Typical values are 
hvnwx = I0rlm e c 2 and hv m i„ = 10~ 8 m e c 2 . The distribution of 
Ush j is set to account for relativistic beaming of the radiation 
in the shock frame as well as for a relatively isotropic compo- 
nent in all frames, from US to DS. This is achieved by a loga- 
rithmic separation of fi s i, in the US direction between /i,/, = 
and fish = lj with 1 -max(^ s /,) < T~ 2 . The fj, S f, < directions 
are chosen as the zeros of a Legendre polynomial, the same as 
the common Gaussian quadrature. A typical division is shown 
in fig. [3] We note that in order to account correctly for the 
relativistic beaming using Gaussian quadrature, for instance, 
would require a much larger number of azimuthal directions 
for high values of T„. The convergence of t he so lutions with 
respect to the resolution is demonstrated in § 15.41 



4.3. Test problems 

The numerical scheme and its implementation were tested 
thoroughly to ensure the results are valid. The tests verified 
a correct description of the different radiation mechanisms 
in steady state problems including, e.g., Compton scattering 
with pair production, bremsstrahlung emission with self ab- 
sorption. We present here only two of the tests, demonstrat- 
ing the suitability of the numerical scheme for dealing with 
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repeated Compton scattering and pair production and annihi- 
lation. The test results are compared with analytic solutions 
or Monte Carlo simulations and are shown to reproduce them 
well. 

4.3. 1 . Comptonization in a cloud of low and medium optical 

depths 

A thin, stationary planar layer of plasma with Thomson op- 
tical depth tj in the z direction (perpendicular to the symme- 
try plane) and a given temperature T is irradiated at one end, 
r* = 0, by a 5 function in v, directed along the z axis, 



/„(/i > 0,t = 0) =I 6(is-v )S(fi- 1). 



(83) 



At r* = tj a free boundary condition, /(/i < 0,r* = tt) = 0, 
is applied. In order to reach the steady state solution for the 
radiation field, the iteration scheme of the radiative transfer 
equation is used until the radiation field converges. 

The results of these calculations are compared with an inde- 
pendent Monte Carlo simulation of the setups usi ng the same 
approximate Compton kernel, as described in § 13.2.11 The 
specific photon flux escaping through the free boundary at 

T* = T T , 



—fidfi, 



(84) 



was calculated for two cases with t = 1, one with tj = 1, 
z>o = 10~ 8 and the other with tj =0.01, vq= 10 -4 . The resulting 
spectra are shown in Figs. |4]and|5]for t t = 1 and t t = 0.01 re- 
spectively. In each figure the results of the code (blue pluses) 
and the Monte Carlo simulation (black lines) are shown. As 
can be seen, there is an excellent agreement between the two 
independent methods for calculating the spectrum of escaping 
photons. 

4.3.2. Pair quasi equilibrium for given T 

This test checks the numerical description of the (integral) 
pair production and annihilation. We use a setup with a given 
Wien spectrum of the radiation field, 



Ip(fl) oc v 



-0/T 



(85) 



For a given f, we find the equilibrium value of x+ = n+/n p for 
which the positron production and annihilation rates cancel 
each other analytically and numerically. A comparison be- 
tween the two values obtained is given in tableQ]for different 
temperatures. 
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FIG. 4. — The specific photon flux leaving a cloud of plasma with f = 1 
and width tt = 1. The radiation entering the cloud has a single frequency 
uq = 10 . The results of the code are marked with blue pluses and the Monte 
Carlo results are shown as a black line. 




FIG. 5. — The specific photon flux leaving a cloud of plasma with f = 1 
and width tt = 0.01. The radiation entering the cloud has a single frequency 
uq = 10 . The results of the code are marked with blue pluses and the Monte 
Carlo results are shown as a black line. 
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TABLE 1 

Equilibrium values of x+, balancing the pair production and 

annihilation rates at different temperatures. 



Note that x + does not necessarily grow with f, since we use 
different densities n for convenience. We obtain an accuracy 
of a few % except for very low temperature, where higher 
resolution is needed in order to account for the exponential 
cutoff near v = 1. The resolution used here is v n+ \/v n = 1.4, 
% = 12. 

5. NUMERICAL RESULTS 

In this section we present the numerical results, solving 
equations (I28li-(l32l self consistently for different values of the 
upstream Lorentz factor T u . We divide the presentation of the 



results into 2 pa rts: The structure (§ I5.lt and the radiation 
spectrum (§ 15. 2t . The structure is the spatial distribution of 
integral parameters such as temperature, velocity (or Lorentz 
factor), pair density and radiation pressure. The spectrum is 
the distribution of radiation intensity at different angles and 
photon energies (at given locations across the shock), mea- 
sured in a specific reference frame. Two important frames 
of reference are the shock frame, in which the solution is a 
steady state solution, and the local rest frame of the plasma, 
which is useful for understanding the interaction between the 
radiation and the plasma. 

5.1. Structure 

The values of TB, T and x+ for T„ = 6, 10, 20 and 30 are 
shown in figures [6] to [TT] as functions of the Thomson optical 
depth for upstream going photons r* [defined in Eq. fiTt l or 
T*/r„. Figures zoomed on the DS region (r* > 0) are sep- 
arately given. The results are calculated for n u = 10 15 cm -3 , 
over regimes where bremsstrahlung absor ption is negligible 
(i.e. they are in the low density limit, see § l2.2.4t . 

The shock profiles can be divided to 4 regions: 

1 . Far upstream - The velocity is constant, while the radi- 
ation intensity and positron fraction grow exponentially 
until they hold a significant fraction of the energy and 
momentum of the flow. 

2. The velocity transition - Here the flow decelerates con- 
siderably, reaching a velocity close to the downstream 
velocity. For RRMS this regime is bound by a sub- 
shock. 

3. Immediate downstream - In the first j3 d l optical depths 
behind the velocity transition the flow approximately 
stays at constant velocity, while the plasma and radi- 
ation are in CE. A gradual cooling by bremsstrahlung 
emission and inverse Compton scattering takes place. 
This region produces the radiation that diffuses up- 
stream and decelerates the incoming plasma. 

4. Far downstream - Further than approximately fi~ d x opti- 
cal depths into the downstream, from where most pho- 
tons can not diffuse upstream. From this point on, a 
slow thermalization takes place accompanied by a slow 
decline in the plasma temperature and photon energies, 
ending when the temperature reaches the downstream 
temperature. The decline in temperature leads first to a 
decrease in positron number, until the pair density be- 
comes negligible compared to that of the original elec- 
trons (x+ < 1) at T ~ 50keV. Then the thermalization 
continues until bremsstrahlung absorption takes over 
and thermal radiation at equilibrium is established. 

We do not solve the equations in the fourth region since the 
solution there is straightforward (the radiation is isotropic and 
in equilibrium with the plasma). Also, note that since the far 
downstream is supersonic, a second sonic point is expected 
in RRMS. This, however, is a stable point with no special 
physical significance. 

Figures [6] and [7] show, for different values of r„, the struc- 
ture of the relativistic velocity T/3 across the shock. It can be 
seen that the deceleration length in units of t* is much larger 
than unity and grows with T u in a manner faster than linear. 
A subshock is obtained at the sonic point, with a discontin- 
uous deceleration of 8(T/3) ~ 0.1. Behind the subshock, the 
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FIG. 6. — The relativistic velocity of the flow F/3 vs. t* /T„ for different values 
of r„, from the US to the subshock (r* = 0). 
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FIG. 7. — The relativistic velocity of the flow V0 vs. t* for different values 
of r„, around the subshock (t, = 0). Notice that the last mean free path on the 
right hand side is influenced by the boundary conditions, however the flow near 
the subshock is not affected by this boundary condition. 
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FIG. 8. — The normalized temperature f vs. T„/r„ for different values of V u , 
from the US to the subshock. 



FIG. 9. — The normalized temperature f vs. r* for different values of F u , 
around the subshock (Notice that the last mean free path on the righ hand side is 
influenced by the boundary conditions). 



velocity approaches its far DS value in a few Thomson opti- 
cal depths. The last optical depth is affected by the boundary 
conditions im posed on the right hand side. This effect will be 
discussed in § 15.4. 21 

Figures [8] and [9] show, for different values of T„, the struc- 
ture of the temperature T across the shock. The far US shows 
an exponential growth of T as a function of t». The tem- 
perature then saturates at a maximum which is approximately 
linear in r„, and then decreases towards the subshock. Be- 
hind the subshock the temperature jumps, reaching a value 
of Tj ump ~ 0.5, which grows with T,,, and then cools with a 
typical distance of a few Thomson optical depths (r*). 

Figures [10] and QT| show , for different values of T u , the 
structure of the positron to proton number ratio, x+, across the 
shock. The growth of x+ as a function of r* when approach- 
ing the subshock is super exponential, and its value reaches 
a maximum a few optical depths behind the subshock. The 
maximal value is approximately linear in T„ (see figure [T6l>. 
Figure [12] shows x+T across the shock, which represents the 
pressure of the positrons and their relative importance in set- 
ting the speed of sound in the plasma, compared to the pro- 
tons. The value of x+T goes above a few hundreds at the sub- 



shock for r„ > 6. 

Figure[T3]shows the ratio of thermal energy flux carried by 
electrons and positrons to the radiation energy flux, F^, vs. 
Tf3/(r u p u ). The energy flux ("taken" from the protons) is 
dominated by thermal and rest mass energy flux of the elec- 
trons and positrons during most of the transition rather than 
by radiation energy flux. The energy is transferred to the radi- 
ation when the flow approaches the DS velocity, and the two 
fluxes are comparable around the subshock. Comparing the 
results at a fixed point (e.g. T = T„/2), this ratio grows with 

r„. 

Figures [14] and Q3] show the relativistic velocity Tj3, the 
temperature T and x+ as a function of the scaled distance 
Zsh = F u n u aTZsh, for r„ = 10. These figures illustrate that the 
shock width is comparable to the upstream Thomson mean 
free path, as z s h is approximately measured in these units. 

5.2. Spectrum 

Figures [17] to [26] show the radiation spectrum at different 
points along the shock profile for the cases r„ = 10 and T u = 
30. The normalization of the intensity and frequency is given 
in Eq. d4TT i. The points of interest are: 
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FIG. 10. — The positron to proton ratio x+ vs. T*/T u for different values of V u , 
from the US to the subshock. 
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FIG. 1 1 . — The positron to proton ratio x+ vs. r* for different values of F„, 
around the subshock. 
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FIG. 12. — The rest frame normalized positron pressure x+f vs. t*/T u for 
different values of r„, from the US to the subshock. 
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1. The upstream - where T = 0.99IY At this point we 
show the spectrum in the rest frame of the plasma (Figs. 
[17]and[T8]for T„ = 10 and r„ = 30, respectively). 

2. The transition - where T = T u /2. At this point we show 
the spectrum in the rest frame of the plasma (Figs. [19] 
and|20]for T u = 10 and T u = 30, respectively), and in the 
shock frame (Figs. Ej] and [22] for r„ = 10 and T„ = 30, 
respectively). 

3. The immediate DS - One Thomson optical depth (r» = 
1) downstream of the subshock. At this point we show 
the spectrum in the shock frame (Figs. |25l and |26] for 
r„ = 10 and r„ = 30, respectively). 

We now give a short description of the main characteristics 
of the spectrum at different locations across the shock. An 
extensive analysis and an analytic description of the results is 
given in section § [6] 

• Upstream: The rest frame spectrum (figs. IT7l and IT8l> 
is strongly dominated by a photon component beamed 
in the US direction, with a typical energy of ~ 3T u m e c 2 , 
and a much weaker, isotropic component with energy 
~ T u m e c 2 . In the shock frame (not shown here), the 
dominant component is beamed in the DS direction, 



with characteristic energy ~ Tlm e c 2 . There is also a 
weaker and not strongly beamed US going component 
with energy somewhat higher than m e c 2 . 

• Transition region: The radiation in this region is ex- 
tremely anisotropic in both the shock frame and the 
rest frame of the plasma. In the rest frame (figs. [191 
l20l the radiation is dominated by a high energy com- 
ponent beamed in the US direction, with a typical en- 
ergy of hv ss Tm e c 2 , where T is the local Lorentz fac- 
tor. An isotropic component, which is much weaker 
in intensity and with typical photon energy similar to 
the beamed component, also exists. In the shock frame 
(figs. [2U l22l the spectrum is composed of a domi- 
nant narrowly beamed component in the DS direction 
with typical photon energy hv ~ T 2 m e c 2 , and of a much 
weaker intensity of US going photons with typical en- 
ergy of hv ~ m e c 2 . 

The spectrum in both frames contains highly beamed 
components. In order to estimate the amount of energy 
carried by the beams, we show in figures [23] and [24] the 
intensity / multiplied by 1 -/i 2 , which for 1 - |/i| <C 1 is 
proportional to the solid angle. In the rest frame, the 
hv re sa Tm e c 2 component dominates the total energy, 
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FIG. 15. — Same as fig. 1141 showing only the DS region. 




FIG. 16. — The maximum value of x+ for different r, ( values. The ap- 
proximation Xf. max = 142r„ is accurate to better than 10% in the range we 
investigated. 

while in the shock frame the energy carried by the US 
going hv s k m m e c 2 photons is comparable to that of DS 
going hv S h ~ T 2 m e c 2 photons. 

• Immediate DS: Figs. |25l and [26] show that the spec- 
trum is composed of two components: a relatively 
isotropic component with hv ~ m e c 2 , and a compo- 
nent narrowly beamed into the DS direction with en- 
ergy reaching hv ^T 2 . 

Figures [27] and [28] show the spectrum integrated over [i 
in the immediate DS. The integrated spectrum is domi- 
nated by photons of energies ~ m e c 2 , but includes a sig- 
nificant high energy tail. The high energy component 
holds 10%-20% o f the total energy flux of the radiation 
and is analyzed in[i 



5.3. Compton scattering and pair production optical depths 

The dominant mechanisms affecting the radiation in the 
transition region are Compton scattering and photon-photon 
pair production. To determine the relative importance of the 
two processes and obtain a handle on some of the important 
physical features of the deceleration mechanism, we examine 
the optical depth for US going and DS going photons in the 



transition region, for the cases r„ = 10 and T u = 30. Figures 
[29l and [30] show the cumulative optical depths for US going 
photons leaving the subshock and reaching the point where 
T = r„/2 as a function of shock frame frequency. It is clear 
that many of the photons with i> s f, > 1 will make it from the 
immediate DS to the middle of the transition, while low en- 
ergy photons v s h <C 1 will be scattered on the way. 

Figures [3J] and [32] show the cumulative optical depths for 
US going photons with v s f, « 1 leaving the subshock, vs. the 
relativistic velocity Tj3 of the flow in the transition region. 
These photons constitute the majority of the photon flux leav- 
ing the immediate DS in the US direction. It can be seen that 
most of the shock profile, up to V ~ 0.9T,,, has a total optical 
of ~ 5 for these photons, most of it due to Compton scatter- 
ing, and order unity optical depth due to photon-photon pair 
production. 

Figures [33] and [34] show the cumulative optical depths for 
DS going photons, starting from the point Y = 3 in the tran- 
sition and reaching the subshock, as a function of v^. Com- 
paring the results for T u = 30 and T,, = 10 we find that the 
optical depth due to both scattering and photon-photon pair 
production are very similar for both values of T u , suggesting 
a common structure and a common upstream going photon 
spectrum in this region. 

Figures [351 and [36l show the cumulative optical depths for 
DS going photons, starting from the point Y = Y„/2 in the 
transition and reaching the subshock, as a function of v S h. As 
was shown earlier, the shock frame radiation in the transition 
region is dominated by photons with energy ~ T 2 m e c prop- 
agating towards the DS. The figures illustrate that the optical 
depth for these photons to reach the immediate DS is less than 
unity. On the other hand, photons with energies around the 
pair production threshold in the shock frame, 0.1 < z> s /, < 10, 
will suffer a strong attenuation due to pair production. 

5.4. Numerical convergence 
5.4.1. Resolution 

The solutio n of t he equations is obtained using iterations, as 
described in § 14.11 Iterations are continued until the changes 
in integral quantities (T, T, x+, P ra d etc.) are less than ~ 1% 
between successive iterations. The resolution used for the so- 
lutions presented in the preceding sub-sections is given in ta- 
ble 

We found that the solutions are modified by ~ 1 % when the 
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FIG. 18. — The plasma rest frame radiation spectrum vie vs - &, for T„ = 30 in 
the US (T = 29.7). 
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FIG. 19. — The plasma rest frame radiation spectrum vlp vs. v, for Y u = 10 in 
the middle of the transition (r = 5). 
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FIG. 20. — The plasma rest frame radiation spectrum £■/=> vs. v, for r„ = 30 in 
the middle of the transition (r = 15). 
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TABLE 2 

The resolution used for the calculations of the shock 

profiles. 
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resolution in t, is increased from St, = 0.2 to St, =0.1, and 
therefore concluded that solutions obtained with either reso- 
lution are satisfactory. The convergence of the solutions with 
respect to the resolution in v^ and /i s /, was tested using solu- 
tions with lower and higher resolutions for r„ = 10. We used 
several properties of the solution to quantify the convergence. 
The solution properties we checked were: 
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shock; 



the temperature immediately behind the sub- 
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10 



I 10" 



10 



10 









H sh =-0.986 






\i h =-0.827 




,- \ 


\i .=-0.108 
\x ,=0.762 

\i ,=0.999 

r sh 




t \ 

: / / « 


H . =1-4.3e- 


05 


'1 V / ' ' ' 




• 


7 V ' ' > 

J u\ \ 11 

/ 'T '** \ n 



10 



10' 



10" 2 

hv Vm c 

sh e 



10' 



10" 



FIG. 24. — The solid angle weighted shock frame radiation spectrum, 
%hlsh,0 /,(1 — A^> f° r ^" = 30 witihn the transition region (r = 15). 




FIG. 25. — The shock frame radiation spectrum £,/,/,/,. => v/? vs. £>,;,, for Y u = 10 
in the immediate DS (t* = 1). 
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FIG. 26. — The shock frame radiation spectrum VstJshfi h vs - ^sh> f° r T» = 30 
in the immediate DS (r» = 1). 
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FIG. 27. — The integrated spectrum in the shock frame, z>,/, J i s h.v, h (Psh)d^i s h, 
at different depths inside the immediate DS (r* > 0) for T,, = 30. 



FIG. 28. — The integrated spectrum in the shock frame, u,i, J/ s ;, => ft (/i,/,)fifyt,/,, 
in the immediate DS (t* = 3) for T„ = 10, T„ = 20 and T„ = 30. 



• The maximal x+ value; 

• Pshjump - the value of the radiation pressure in the shock 
frame at the subshock; 

• — T*(r/? = 5), the normalized optical depth upstream 
of the subshock at which the Lorentz factor drops by 
-half; 



• -T*, n i - the normalized optical depth upstream of the 
subshock at which the US ev olutio n becomes nonlinear 
(see detailed explanation in § 16.31 ; 

The value of r* .„/ is very sensitive to small changes in resolu- 
tion, since it is set by the exponential decay of the number of 
photons arriving from the immediate DS. However, its exact 
value does not affect significantly the structure of the decel- 
eration region. We use it here merely as a stringent test of 
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FIG. 29. — Cumulative optical depth of US going photons from the subshock to 
r = r„/2 vs. shock frame frequency u s /,, due to Compton scattering and photon- 
photon pair production, T,, = 10. 




FIG. 30. — Cumulative optical depth of US going photons from the subshock to 
r = r„/2 vs. shock frame frequency £>,/,, due to Compton scattering and photon- 
photon pair production, T,, = 30. 




FIG . 31 . — Cumulative optical depth of US going photons with i> s i, = 1 leaving 
the subshock vs. T/3/(T u /3u)t due to Compton scattering and photon-photon pair 
production, Y u = 10. 
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FIG. 32. — Cumulative optical depth of US going photons with i> s t, = 1.1 leav- 
ing the subshock vs. r/3/(r„/3„), due to Compton scattering and photon-photon 
pair production, Y u = 30. 



numerical convergence. 

The changes in the values of the test parameters as a func- 
tion of resolution are given in fig. [37] The results were ob- 
tained using lower and higher resolutions in v and /i, and are 
presented as a function of N v x N^, the product of the number 
of discrete values chosen for fi and for v. The numerical er- 
ror around the nominal resolution used in our calculations is 
few percent, except for the most sensitive parameter, T* >n i, for 
which the numerical error is around 10%. 

5.4.2. Changes in the length of the DS 

In order to verify that the boundary conditions imposed on 
the DS edge of the shock do not have a significant effect on 
the final results, around the subshock and in the shock transi- 
tion region, we compare the results shown above to the results 
obtained with a solution including a longer DS region behind 
the subshock. We are limited in extending the DS because 
of numerical problems, caused by the proximity to a second 
sonic point. For this reason we extend only the DS of the cal- 
culation for r„ = 30, from t* = 7 in the calculations presented 
above to r* = 10. The changes in integral quantities resulting 
from this modification of the DS region length are of order 
of a percent. The only quantity which changes by a larger 
amount, ~ 3%, is t*.„/. The temperature and velocity profiles 



obtained in the two calculations are compared in figs. [38land 



6. SIMPLIFIED ANALYTIC MODELLING OF RRMS 
STRUCTURE 

The key to understanding the structure of RRMS lies in the 
understanding of the behavior in the immediate DS. In our 
qualitative analysis of the immediate DS of RRMS, § 12.3.31 
we have argued that the immediate DS photon-electron- 
positron plasma should be close to Compton pair equilibrium 
(CPE). This enabled us to demonstrate that the temperature 
in the immediate DS is expected to be T s ~ QAm e c 2 , and that 
the immediate DS should be sub-sonic. These results are con- 
sistent with the numerical resul ts presented in § [5] We first 
discuss in some detail in § 16.11 the accuracy of the CPE ap- 
proximation for the description of the immediate DS. 

Once the immediate DS is understood, a simple estimate 
of the spectrum of photons emanating from this region in the 
US direction leads to an understanding of the transition (de- 
celeration) re gion, and of th e asymptotic (far) US. These are 
discussed in § !6.2l and § 16.31 

The flow downstream of the immediate DS is smooth and 
NR. As the plasma flows away from the shock transition, it 
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FIG. 33. — Cumulative optical depth of DS going photons from the point T = 
3 to the subshock vs. u s /,, due to Compton scattering and photon-photon pair 
production, T,, = 10. 




FIG. 34. — Cumulative optical depth of DS going photons from the point T = 
3 to the subshock vs. *>,/,, due to Compton scattering and photon-photon pair 
production, r„ = 30. 
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FIG. 35. — Cumulative optical depth of DS going photons from the point V = 
r„/2 to the subshock vs. £> s i,, due to Compton scattering and photon-photon pair 
production, T u = 10. 
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FIG. 36. — Cumulative optical depth of DS going photons from the point T = 
r„/2 to the subshock vs. £>,/,, due to Compton scattering and photon-photon pair 
production, r„ = 30. 
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FIG. 37. — A summary of numerical convergence tests. Results obtained 
with different resolutions in v and fi are shown for T M = 10, as function of 
N v X Nfi, the product of the number of discrete values chosen for fi and for v. 
The y axis shows the values of Tj„ mp , x+jnax, PshJump* ~ T * (T/8 = 5) and — r*.„/ 
(see text for definitions), divided by their values obtained using the resolution 
given in table[2] The different sets of {N^Nv} shown are: {8,26}, {8,41}, 
{13, 16}, {20,41} and {13,60}, while the reference is {13,41}. 



slowly produces the photon density needed for thermal equi- 



librium, eventually reaching the asymptotic DS thermal equi- 
librium conditions. This "therma lizat ion" phase is discussed 
in § 16.41 Finally, we discuss in § 16.51 the high energy photon 
"beam" propagating from the transition region into t he D S, 
and comment on the behavior in the r„ — > oo limit in § 16.61 

6.1. Immediate DS 

Let us examine the accuracy of the CPE approximation, re- 
lating the temperature, the number of positrons and the spec- 
trum of the photons. In the top panel of figure [40] we com- 
pare the average photon energy (i>) e ff and the «;/n 7 , e ff ratio 
obtained in the immediate DS of the T u = 20 solution, with 
those expected at CPE, (i>) e s = 37\ and 



»7,eff 



dxx e 



2 e -V* 2 +f- 2 _ K ^ T ') 



T- 



(86) 



where K2 is the order 2 second kind modified Bessel function. 
The effective number of photons and the average energy per 
photon were calculated using the spectrum around the maxi- 
mum of I u in the rest frame of the plasma, i/ pea k. Specifically, 
the numerical values shown for n 1>e n and (i>) e ff are the num- 
ber of photons in the energy range [f pea k/10, 10^ pea k] and 
their average energy, respectively. The values of (£) e ff and 
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FIG. 38. — f vs. r* for T,, = 30, obatined using a DS optical depth of r* = 7 
(dashed line) and r* = 10 (solid line). 



FIG. 39. — T/3 vs. r* for r u = 30, , obatined using aDS optical depth of t„ = 7 
(dashed line) and r» = 10 (solid line). 



of n//n 7 . e ff are shown as functions of f in the vicinity of the 
sub-shock. Immediately downstream of the sub-shock these 
values are far from those expected for CPE, and they approach 
the CPE values away from the sub-shock. The figures show 
a systematic deviation from CPE. This is expected, sinc e the 
high energy DS photon beam, discussed in detail in 16.51 car- 
ries a significant fraction of the energy and is very weakly 
coupled to the plasma, due to suppression of the cross sec- 
tions at high photon energy. 

The lower panels of fig.[40]compare the values of V, P m d,sh, 
and x+ obtained in the numerical solution, with those obtained 
under the CPE approximation (note, that under the CPE ap- 
proximation the conservation eqs., eqs. d28l . (129k and (l30i l. 
allow one to determine T, P ra d,sh, and x+ as a function of T). 
Here too, the solution deviates from the CPE predictions im- 
mediately downstream of the sub-shock, and approaches the 
CPE prediction away from it. We conclude that the CPE ap- 
proximation yields estimates of the global flow variables (T, 
Pmd,sh, and x+) which are accurate in the immediate DS to 
within tens of percent. 

6.2. The transition region 

The transition or deceleration region is the region in which 
the energy and momentum flux of the US plasma are trans- 
ferred to the radiation and to the e + e~ pairs. The behavior in 
the r„ ^> T 3> 1 regime may be understood using the follow- 
ing arguments. 

1 . The photons decelerating the plasma originate in the 
immediate DS and have a shock frame energy of~ m e c 2 
and a rest frame energy of ~ Tm e c 2 . Since the imme- 
diate DS temperature is T s ~ 0Am e c 2 (see § 16. U . the 
characteristic shock frame energy of these photons is 
hv ~ 3T S ~ m e c 2 . Across the transition region, these 
photons dominate the energy density in the rest frame 
of the plasma, where their energy is ~ Tm e c 2 (see e.g. 
figure |23T>. 

2. The upstream going photons decelerate the plasma by 
Compton scattering, and by pair production interac- 
tions with photons, that are generated within the tran- 
sition region either by Bremsstrahlung emission or by 
inverse Compton scattering (upstream going photons 
that are back-scattered by the downstream flow). The 



three processes similarly contribute to the deceleration, 
as explained in point[8]below. 

T ~ Tm e c 2 . Both Compton scattering and photon- 
photon pair production generate electrons/positrons 
with characteristic energy ~ Tm e c 2 , driving the plasma 
temperature to ~ Tm e c 2 . 

Pairs produced in the deceleration region drift with the 
plasma all the way to T ~ 1 without annihilating, due to 
the high temperatures that reduce the annihilation cross 
section (ex \og2T /T) and to the ex T~ 2 suppression of 
the collision rate. 

The plasma rest frame energy density is dominated 
by pairs rather than protons (n e + n+)(m e c 2 + 3T) > 
n p m p c 2 . Once V <C T„, most of the energy flux is car- 
ried by radiation and pairs. The pairs carry a significant 
fraction of the energy flux (see figure [13), hence their 
energy in the rest frame dominates over the protons rest 
mass. 



6. A significant deceleration of a fluid element, T — > T/2, 
requires that the number of Compton / pair production 
interactions occurring within it be similar to the num- 
ber of leptons within it. A change of factor 2 in V cor- 
responds, in the plasma rest frame, to an acceleration 
to velocity j3' = 0.6 towards the US. This requires a 
momentum transfer of ~ TT jc to each lepton (recall, 
that for r <C T,, the plasma energy density is dominated 
by pairs). This is similar to the momentum transfer by 
Compton scattering or pair production interaction of a 
typical US going photon, for a plasma rest frame tem- 
perature of ~ 2Tm e c 2 . 

1 . The optical depth for typical US going photons between 
r — > T/2 is At ~ 1. The number flux of US going 
photons is similar to the sum of number fluxes of typi- 
cal DS going photons and pairs (pairs are downstream 
going). The similarity between the number densities 
and the fact that DS going pairs undergo ~ 1 interac- 
tion between T —> T/2 implies that US going photons 
roughly interact once as well. 

8. The Thompson optical depth in the range T — > T/2, 
is roughly At, ~ V 2 , wether the deceleration is 



21 




0.16r 



0.14 



0.1 - 



0.08- 



0.06 



0.2 



0.3 



0.8 
0.7 

0.6 

_ °- 5 

c^0.4 

c 
0.3 

0.2 

0.1 



1.1 0.2 0.3 0.4 0.5 0.6 
T/m c 2 



0.7 0.8 



4000 1 

3500 
3000- 
2500 
< + 2000 
1500 
1000- i 
500 / 

o L 



0.1 0.2 0.3 0.4 05 0.6 
T/m <? 



0.7 0.8 



0.4 0.5 

T/m cf 



0.6 



0.7 



FIG. 40. — Comparison of the numerical solution in the immediate downstream with results obtained assuming CPE, for the case T u = 20. The blue dash-dotted 
curves show the values obtained assuming CPE, while the green dashed and red solid lines show the numerical solution behind the sub shock. Green dashed lines 
show the results in the first Thomson mean free path behind the sub shock and the red solid lines show the results deeper into the downstream. 



due to Compton scattering or due to pair creation 
on Bremsstrahlung generated photons. This implies 

r»(r)~r 2 . 

In the range Y —$ Y/2 there is ~ 1 interaction per lep- 
ton crossing (see point |6j, The Thomson optical depth 
required for a single Compton scattering is 8r* tS cM ~ 
IT ~ T 2 due to the KN correction to the cross sec- 
tion. Similarly, the Thomson optical depth required for 
a single pair production on a "returning" (downstream 
scattered) photon is ST* >re t ~ r 2 n 7rer /n/ ~ r 2 , where 



■ «/ is the number density of returning photons. 



Bremsstrahlung generated photons with energy hv ~ 
m e c 2 /T have a large optical depth for pair creation on 
the US going typical photons, since they do not suffer 
a suppression to the cross section. The number of these 



photons, produced up to a given point in the transition 
region is given by 



n lt ff = Qjjf 



ksh 



(87) 



where Q-,.// is the production rate of photons that are 
able to upscatter to hv ~ m e c 2 /r in the rest frame of the 
plasma and &,/, is the shock frame distance over which 
r changes significantly. The Thomson optical depth re- 
quired for producing enough photons to decelerate the 
plasma, n 7 ,// ~ «/, is thus 



St, 



'*,// ' 



(88) 



a e g(T,i>)A U s 
where 10 < g(T, v) < 20 is the the Gaunt factor at high 
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temperatures and low frequencies hv/(m e c 2 ) ~ T~ l and 
Aj/5 ~ 5 is a logarithmic correction accounting for pho- 
tons that are produced at low energies and upscatter 
to the required energy by the available ~ 1 number of 
Compton scatterings on the thermal electrons. 

Two conclusions can be drawn. First, the Thomson 
cross section needed to decelerate a Lorentz factor V 
is r* ~ r 2 . Second, all three processes discussed in this 



point are comparable. Simply taking V 



(where 



r* is measured from the subshock) results in a qualita- 
tively good fit to the numerical results, as can be seen 
in fig. [41] It is evident that the deceleration, when ap- 
proaching the subshock, has a universal structure for 
different T u values. 

The following additional properties are implied by the 
above considerations. 

• x + ~ (r„/r 2 ) x (m p /m e )/8 when T -C T„. This follows 
from conservation of momentum flux, and the signifi- 
cance of pairs in the flux. 

• Most of the photons resulting from Compton scattering 
will propagate to the immediate DS without undergoing 
further interactions. The optical depth for scattering of 
photons originating from scattering into the DS direc- 
tion is negligible since they have shock frame energy of 
~ T 2 m e c 2 and suffer a ~ T~ 2 attenuation in interaction 
rate. The optical depth for pair production is of order 
unity. This can be seen by the fact that the cross sec- 
tion and target photons for pair production are similar 
to the Inverse Compton cross section and target photons 
of the e + and e~ in the deceleration region. In fact, the 
total optical depth is < 1 , as can be seen in figure l36l 




FIG. 41. — The relativistic velocity Tf3 for different values of V u and the 
simplified anlytical result for the structure T/3 ~ v / — t\T (bold line) vs. r* . 



6.3. Far US 

As was shown in 15. 21 the radiation, as seen in the rest frame 
of the far US plasma, is strongly dominated by a beamed 
[fi resl « — 1 + 1/r 2 ], radiation field with photons of typical 
energy of several times T u m e c 2 . To understand the main phys- 
ical properties of this region, it is useful to approximate the ra- 
diation field as a delta function in energy and direction, going 
in the US direction. The asymptotic solution for such a radi- 
ation field can be easily found to be an exponential growth of 




FIG. 42. — T,, = 10 far US exponential growth of integral quantities. The 
Cyan lines s how e XT * and e 2Ar * , which are the simplified model exponen- 
tials (see eg. 1891 expected for P, F and T (A) and x+ (2A). 

the parameters P ra d,sh, Frad,sh, r„/3„-r/3 and T, with the same 
exponent, X as , and an expon ential growth of x+ with an expo- 
nent 2X as fsee lSagivl (120061) 1. An approximate value for \ as 
is given by 



K 



*0.28- 



1 



r„/3«n„(Tr 



(89) 



where a c is the total cross section for the photons in the rest 
frame of the plasma. The results of the numerical calculations 
are shown in figure|42]to agree with the expected exponential 
growth for the case T u =10, where the rest frame dominant 
frequency is hv rest f=a 36m e c 2 . 

The solution deviates from exponential growth when the 
temperature approaches m e c 2 as, for example, the Compton 
cross section changes significantly. For r u = 10 the transi- 
tion from linear to non-linear evolution occurs at t* „/ w -320, 
while for Y u = 20 it occurs at t*„/ f=a -890. r* t „; grows with 
r„ - 1 (energy per proton) in a manner faster than linear. 

6.4. FarDS 

This region is characterized by an almost constant veloc- 
ity and a slow growth in photon number that lowers the tem- 
perature. It can be divided into two regions: T > 50 keV, 
where x+ > 1 and electron-positron annihilation takes place, 
and T < 50 keV, where the number of positrons is small and 
they play no significant role. 

Let us first consider the x+ > 1 region. The low temperature 
limit of eq. d86l ) yields 



*7,eff 



(r<i): 



= 2, I' 



st 3 / 2 



(90) 



Using arguments similar to those used for the immediate DS 
estimates, we can write an equation for the evolution of pho- 
ton number 



1 dtt 7 ,eff(T*) gf f (f)A eff (f) 



h/(t*) c/t* 



PdVT 



(91) 



and assuming that most of the energy flux is already in the 
radiation we can approximate n^ e ffT = Const. Using this as- 
sumption with Eqs. (|90T > and (|9T| l we obtain 



dT_ =2 hr g ff (T)A eff (T) ^_ uf 
dr* V 8 Pd 



(92) 
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The flow reaches n e + « n p (x+ sa 1) when f sa 0.06, with a 
weak dependence on r„. From Eq. J92l ) we see that the length 
scale is set by the lowest temperatures in this range, and rea- 
sonable parameters yield ~ 10 5 optical depths required for the 
positrons to annihilate. 

In the the x + < 1 region, a gradual increase in photon 
number lowers the temperature until thermal equilibrium is 
reached. The length scale for this process is 



■ PdC 



Qeff.d 



(93) 



where n 7 , eq ~ obbT] /2.8 is the thermal equilibrium photon 
density and 

Qeff.d ~ geff,dh-eff,d n2 d< J T c l vt( is the photon generation rate 
in the DS. An estimate for Lj yields 



(a T n d ) ' 



3 x 10 6 (3&) 



geff,d^eff,d 
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r 3/4 -1/8 
1 u.2"u.l5 ■ 



(94) 



In this region, the temperature drops as a power law T(r*) ex 
(r* - To)" 2 , as was shown in lKatz et al.l d2010l) . 

6.5. The high energy photon component beamed in the DS 

direction 

As can be seen in Figs. [25] and [26] the immediate DS has 
a high energy photon component beamed in the DS dir ectio n. 
We use below the simplified analysis presented in § 16.21 to 
derive the characteristics of the spectrum of this beam. 

The photons in this beam originated from the immedi- 
ate DS, propagated into the transition region and then were 
Compton scattered once before returning to the DS. Photons 
that were scattered at a point with Lorentz factor Y return 
to the DS with an energy boosted to ~ Y 2 m e c 2 and within 
a beaming angles 9 ~ T" 1 . 

Denote the shock frame intensity of US go ing p hotons with 
typical energies /o(— r*). Conclusion [7] in § 16.21 leads to the 
equation 

/o(D = er/o(r/2) (95) 

where er ~ 1/3 is related to the exact total optical depth for 
typical US going photons from Y/2 to Y. Assuming that the 
fraction of photons that scatter is constant with Y, the resulting 
intensity emitted at Y, Ib(Y) will be 



/ fl (r)w4e r /i,(r/2) ) 



(96) 



since the photons gain a factor of ~ Y 2 to their energy when 
scattered at Y. The scattered photons are beamed into a cone 
with an opening angle T -1 in the DS direction. Since the 
losses of the scattered photons on the way to the immediate 
DS are less than a factor of 2 and depend weakly on the an- 
gle and energy of the photon, we find that the spectrum of the 
high energy beam Ib can be approximately described as 



v S hh{i> S h,0 S h) oc v°? ©{9,1 - vl )Q(v max -v s h), 



(97) 



where v max sw Y 2 and oc\ sw log 2 (4er)/2 is close to zero, and is 
equal to zero when 4er = 1 . 

We next verify that this analysis complies with the numer- 
ical results (results shown for Y u = 20 calculation). Fig. [43] 
shows the shock frame intensity of a beam with #,/, « 10~ 2 
with different v s j, along the shock, vs. Y 2 /v^. We see that the 
intensity is mostly contributed by the part in the flow in which 
T 2 fa 200P s h, as the physical picture requires. Fig. l44l shows 



the shock frame intensity immediately after the subshock, at 

1/2 

different v S h, as a function of v S h v Jh ■ We see that the struc- 
ture of the beam is such that the different energies are beamed 
according to Eq. d97b . 

6.6. The Y u — > oo limit 

Based on the results for r„ < 30 and the analysis above, it 
appears that for Y u — > oo, T(r*) and Y(t*) approach asymp- 
totic profiles in the regime where Y <C Y u and x+ ^> 1. In 
particular, T ~ m e c 2 in the immediate DS and Y ~ r 2 in the 
transition region. However, we have also seen that the high 
energy beam becomes more dominant as Y u grows. The struc- 
ture of the shock, particularly the immediate DS, may be dif- 
ferent if the high energy beam becomes the dominant carrier 
of momentum and energy of the radiation. Unfortunately, a 
full calculation of very high r„ shocks is beyond our current 
numerical capabilities, and requires further investigation. 

7. NR RMS REVISITED 

In this section we briefly describe a preliminary application 
of the code to NR shocks. Our numerical scheme was de- 
signed and optimized for the solution of the relativistic prob- 
lem, and is not efficient and easy to use for NR problems. 
The main difficulties are 1 . Solving the momentum and en- 
ergy conservation equations for the velocity and temperature 
of the plasma is problematic due to the negligible contribu- 
tions of the thermal energy and pressure. 2. Radiation field 
convergence requires a large number of iterations, roughly 
one iteration per single Compton scattering, implying oc fi~ d 2 
iterations. 

A different scheme for finding the plasma temperature and 
velocity and a different boundary condition in the far DS were 
used for solving the NR problem: 

• The temperature was set to the local CE value calcu- 
lated from the radiation field and the velocity was found 
by solving the momentum conservation only. This ap- 
proximation is justified in the case where n 1 /n e ^> 1, 
where n 7 is the number density of photons, which holds 
in the transition region (when the energy density of the 
radiation is a fair fraction of the flow) and the down- 
stream of a NR RMS. Convergence required that the 
temperature be set to a value that is slightly smaller than 
the actual CE value. 

• The following downstream boundary condition was 
used. The radiation field in the upstream direction 
IpXjish < 0,r» = max(r)), was set to represent the ra- 
diation field at a chosen point in the downstream. This 
was done by assuming a Wien spectrum with a temper- 
ature lower than T s [see Eq. ([T71il. and an intensity that 
satisfies the equilibrium at the DS velocity as expected 
in the DS well behind the velocity transition. 

The radiation transport is solved similarly to the relativistic 
case. 

Figures a preliminary solution for a shock with upstream 
energy per proton e = 50 MeV and a very low density (n„ = 
10 6 cm" 3 ), which ensures that bremsstrahlung absorption re- 
mains unimportant until after the velocity has already reached 
its downstream value. In the calculation shown here, absorp- 
tion is everywhere unimportant, since it does not reach the 
downstream temperature. We stress that the resulting solu- 
tion contains a limited region of optical depth ~ /3j behind 
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the velocity transition, hence the temperature profile may not 
correctly represent the actual solution. 

Figure [45] shows the structure (T/3, T and P) of a shock 
with e = 50 MeV as a function of r*. The dotted black line is 
the analytic solution for Tj3 obtained by Weaver 1976 [equa- 
tion (5.10) there], using gq = 0.56(77- for the average Compton 
cross section, suitable for hv « Q.5m e c 2 typical photon energy 
in the transition region. Weaver's solution deviates from the 
numerical solution near the immediate DS. This is due to the 
lower average photon frequency there, compared to that in the 
transition region, which leads to an increase in ac in the (more 
accurate) numerical calculation. 
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with average Compton cross section ac = 0.56or. 



Examining the radiation spectra obtained in our numerical 
calculations, fig. |46] we are able to verify the validity of two 
of Weaver's assumptions. First, it is clear that the spectrum at 
each point along the shock is dominated by photons within a 
narrow energy range. Second, the anisotropy of the radiation 
is of order (3, which is the expected anisotropy due to diffusion 
of the radiation. 

To conclude, the preliminary solution found using our nu- 
merical scheme is consistent with Weaver's results. In addi- 
tion, the detailed spectra support the validity of Weaver's ap- 
proximations regarding the radiation spectrum. The fact that 
the results for NR shocks are in agreement with previous work 



supports the validity of the numerical scheme. 

8. DISCUSSION 

We have calculated and analyzed the structure of relativistic 
radiation mediated shocks (RRMS). A qualitative discussion 
of the shock physics was presented in § including analytic 
estimates of the deceleration and thermalization length scales 
of non-re lativistic (NR) RMS (equations ( fTTT i, (fl4l >; see fig- 
ure |222] for a schematic shock structure description) and of 
the immediate DS temperatures of both NR RM S [eq. ( fTTb l 
and RRMS [eq. (l26ll. We have also shown (in § l23Jb that 
the immediate DS of RRMS is expected to be subsonic, and 
concluded that the structure of RRMS must include two sonic 
points. 

In section §|3]we derived a dimensionless form of the equa- 
tions describing the conservation and transport equations de- 
termining the structure of the shock, and described in detail 
the radiative processes included in our treatment and the ap- 
proximations we used. In section § [4] we presented a novel 
iteration scheme for numerically solving the equations, and 
demonstrated its validity by applying it to several test cases. 
In section §[5] we have presented numerical solutions for the 
profiles and radiation spectra of RRMS, for upstream Lorentz 
factors r„ in the range of 6 to 30. The main results obtained 
are described below. 

[1] Structure and radiation spectrum. In § !5.1l we showed 
that the structure of RRMS can be divided into four regions, 
from upstream (US) to downstream (DS): The far US, the 
transition region, the immediate DS and the far DS. The far 
US is characterized by a velocity close to the US velocity and 
a radiation energy-momentum flux much smaller than that of 
the US plasma. The transition region is where the velocity 
(T/3) changes significantly, approaching Tf3 ~ 1, while the 
momentum and energy fluxes are transferred to the e + e~ pairs 
and to the ra diati on. In both regions, the radiation spectrum 
(shown in § I5.21 > is dominated in the plasma rest frame by 
US going photons with energy of a few times Tm e c 2 . In the 
shock frame the radiation is dominated by DS going photons, 
beamed into a cone with opening angle ~ T -1 , and a typical 
energy T 2 m e c 2 . In the far US the temperature grows expo- 
nentially with t* towards the downstream (fig. l42l . until it 
reaches ~ m e c 2 (r* is the Thomson optical depth for photons 
moving towards the upstream). The temperature then contin- 
ues to grow at a slower rate until it reaches T /(m e c 2 ) ~ T in 
the transition region, and then decreases, approximately fol- 



25 




10" 



10 



10" 



10" 











(^=-0.932 

..H sh =-0.239 
,.,.,M. sh =0.239 






^=0.932 


/ V\ 

/ V 






I 1 




\l 




V 



hv/m c 



10" 



10 2 

hv/m c 



10 



10" 











0° 




H =-0.932 
|i h =-0.239 




o- 2 




u h =0.239 
u h =0.932 




o- 4 






\ 


o- 6 




\ 



10"' 



hv/m c 



10 



FIG. 46. — Spectra of the radiation in the shock frame along the shock profile for e = 50MeV. Upper left: far upstream (/3 = 0.99/3,,), upper right: inside the 
velocity transition (/3 = 0.5/3,,) and lower: In the immediate downstream (t, = 37). 



26 



lowing the deceleration, T/(m e c 2 ) ~ Y (fig. IS). 

The transition region ends at a subshock, possibly mediated 
by plasma instabilities, with a velocity jump of S(Yf3) ^0.1 
and a slight increase in temperature, to 0.4 < T /(m e c 2 ) < 0.9 
for 6 < r„ < 30 (see Figs. and 01. The immediate DS, fol- 
lowing the subshock, is characterized by a small change of 
velocity, approaching the DS value within ~ 2 Thomson op- 
tical depths, and a temperature that decreases on a scale of 
a few Thomson optical depths to T /(m e c 2 ) ~ 0.25. The ra- 
tio of positron density to proton density in the immediate DS 
reaches a maximum of ~ 140r„ (see Fig. [TTJ, approximately 
when the temperature crosses T /(m e c 2 ) ~ 0.3, and then de- 
creases. The radiation spectrum in the immediate DS is dom- 
inated by a relatively isotropic component with hv ~ 37\ but 
a fraction of 1 0% - 20% of the energy flux is carried by a high 
energy photon tail, strongly beamed towards the DS, with a 
cutoff at ~ Y 2 m e c 2 and a nearly flat spectrum, vF v oc v° (see 
Figs.|26]to|2g. 

[2] Optical depth s du e to Compton scattering and pair 
production: In § 15.31 we showed that the optical depth of 
typical photons (hv ~ m e c 2 ) leaving the subshock in the US 
direction is a few. The optical depth is provided by both 
Compton scattering and pair production, the latter having a 
somewhat smaller contribution (see Figs. l3Tlandl32l. Photons 
with much smaller energies are scattered close to the imme- 
diate DS and do not reach the transition region (see Figs. |29l 
and 13 Oil. Typical DS going photons from the transition region 
(with shock frame energy ~ T 2 m e c 2 ) undergo very few inter- 
actions on the way to the immediate DS (see Figs. [33ltol36ll. 
[3] The importance of e + e~ pairs. In figure Qj] we show that 
the pairs produced along the shock transition and in the imme- 
diate DS play an important role in decelerating the US plasma. 
The energy flux removed from the protons is dominated by 
pairs over radiation during most of the transition, and the ratio 
between the two becomes larger as T u grows. The pair energy 
flux is dominated by thermal energy flux since the transition 
region temperatures are relativistic (T > m e c 2 ). 

We find several characteristics of the structure of RRMS, 
which are qualitatively different from those of NR RMS. 

1 . The Thomson optical depth of the transition region is much 
larger than unity, is dominated by pairs, and grows with T u in 
a manner faster than linear. However, the actual (KN cor- 
rected) optical depth (including pair production) for a typical 
photon crossing the shock is of order of a few. 

2. The temperatures of the pair plasma within the transition 
region are relativistic, T > m e c 2 . 

3. The relativistic shock structure includes a sonic point cross- 
ing, in which the flow changes from supersonic to subsonic. 
We find that this sonic point must be a sub-shock mediated by 
processes not included in our calculation, which operate on 
a scale much shorter than the radiation mean free path [e.g. 
plasma instabilities, see eq. ([T9|>1. 

4. e + e~ pairs carry most of the energy and momentum flux in 
the transition region 

5. In RRMS a fair fraction of the energy density in the imme- 
diate DS is carried by a nonthermal tail of high energy pho- 



tons, where in the DS of NR RMS the radiation is in CE with 
the plasma. 

We developed in §[6] an analytical understanding of the key 
features of the shock structure and radiation spectrum. Sev- 
eral points should be highlighted. 

[1] Immediate DS. The key to understanding the structure 
and radiation spectrum of RRMS is the understanding of the 
immediate DS. The immediate DS of RRMS is close to CPE 
(see Fig. l40l . which, due to the fast increase of the number 
of pairs with temperature, sets the temperature to a large frac- 
tion of m e c 2 (Katz et al. 2010). The large amount of positrons 
and the high temperature imply a relativistic speed of sound 
in matter (3 SS ~ l/v3, and combined with the low velocity in 
this region that quickly approaches its DS value (3d < 1/3, a 
subsonic regime is inevitable. The immediate DS acts as the 
supplier of photons directed towards the US, which deceler- 
ate the incoming plasma through Compton scattering and pair 
production. 

[2] Deceleration region. For T -c T„, we find r(-r») w ^/-t* 
(see fig. |4T1 >. where the subshock is located at r* = 0. This be- 
havior is due mainly to the KN scaling of the cross sections, 
and to the fact that the optical depth for US going photons is 
of order few. This approximation follows closely the numeri- 
cal results up to T w r„/2. 

[3] High energy photon beam. The immediate DS has a 
high energy photon component narrowly beamed in the DS 
direction, with a nearly flat power-law like spectrum, vl v oc v° 
and an energy cutoff at ~ Y 2 m e c 2 . The photons in this beam 
originated from the immediate DS, propagated into the tran- 
sition region and then were Compton scattered once, before 
returning to the DS. Photons that were scattered at a point 
with Lorentz factor Y return to the DS with an energy boosted 
to ~ Y 2 m e c 2 and within a beaming angle 9 ~ Y~ l . An ap- 
proximate description of the resulting spectral and azimuthal 
structure of the beam is given in Eq. d97l i. The total optical 
depth for these photons to reach the immediate DS is small, 
and they carry 10% -20% of the energy flux in the immediate 
DS. The beam is only stopped far into the DS, producing pairs 
on low energy photons. 

[4] Far US. In the far US, P md ^ F radiSh , Y u j3 u -Yj3 and T, all 
grow exponentially with r* with the same exponent, X as given 
in eq. d89l . while x+ grows exponentially with an exponent 
2\ as (see fig. l42l . 

[5] Thermalization length scale. The thermalization length 
is much longer than the shock transition, both in terms of 
Thomson optical depth and in real distance. Thermal equi- 
librium is reached ~ 10 6 Thomson optical depths into the DS. 

Finally we showed for completeness in §|7]the preliminary 
results of a detailed calculation of the structure of a NR RMS 
including full radiation transport. The results are consistent 
with previously published ones, and support the validity of the 
numerical methods we use and of the diffusion approximation 
used for solving the problem in earlier work. 

This research was partially supported by Minerva, ISF and 
AEC grants. 
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APPENDIX 
A. NOTATIONS FREQUENTLY USED IN THE PAPER 



Subscripts, superscripts and miscellanea 

as : Asymptotic upstream behavior 
US : Upstream 
DS : Downstream 
CE : Compton equilibrium 
CPE : Compton-Pair equilibrium 
d : Asymptotic downstream (postshock) value 
dec : Deceleration 
e : Electron value 
p : Proton value 
NR : Non-relativistic 
pi : Plasma value 
rad : Radiation field value 
sh : Shock frame value 
rest : Plasma rest frame value 
u : Asymptotic upstream (preshock) value 
7 : Photon value 
+ : Positron value 
a (hat) : Normalized units 



Symbols 

a\,b = '■ Radiation constant 
a c : Compton scattering cross section 
F ra d ( er g s cm" 2 s"') : Radiation energy flux 
/z(ergs s) : Planck's constant 

I(tt, i/)(ergs cm -2 s _I str -1 Hz" 1 ) : Specific intensity of radiation field 
7/(17, z/) (ergs cm" 3 s" 1 str" 1 Hz" 1 ) : Emissivity coefficient 
l'(cm) : Photon mean free path 

n ( ,,n + ,n,,n 7e ff(cm" 3 ) : Number density of electrons, positrons, ions (protons) and typical photons 
n u (cm -3 ) : Upstream proton (and electron) number density 
P (ergs cm" 3 ) : Pressure 
Prod (ergs cm" 3 ) : Radiation pressure 
2 + (cm" 3 s _1 ) : Net rate of positron production 
T (erg) : Electrons & positron temperature 
T = T /m e c 2 
T Qz , T zz (eigs cm" 3 ) : Components of stress-energy tensor (energy and momentum fluxes, 

respectively) 
x+ = n+/ni : Positron fraction 
z(cm) : Length along flow direction 
X(ri,i/)(cm _1 ) : absorption coefficient 
a e : Fine structure constant 
fi = V 1 - T" 2 : Flow velocity (units of c) 
F : Flow Lorentz factor 
T„ : Upstream flow Lorentz factor 

7 e fh '■ Lorentz factor associated with random motion of e + and e" 
(5=1- Tf3/T u f3 u : Asymptotic deceleration parameter 
e sc , (ergs) : Radiation emission cutoff energy due screening 
£ : Riemann's zeta function 

r\ = exp(~7£) = 0.5616 (where je — 0.5772 is Euler's constant) 
\(.ff) • Correction factor for bremsstrahlung emission 
Xd (cm) : Debye length [= y/T /A-Ke 2 {n e + n+)] 
[i : Cosine of angle relative to positive z-axis (flow direction) 
v (Hz) : Photon frequency 
v = hv/m e c 2 

a c (cm 2 ) : Total Compton scattering cross section 
a 77 (cm 2 ) : Cross section for 77 — ^ e + e~ pair production 
ot (cm 2 ) : Thomson cross section (87rro/3) 
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t* : Thomson optical depth for upstream-going photons, given by r* = [T(l + f3)(n e + nf)ajdz s h. 
r„ = at the subshock and grows towards the downstream. 

B. COMPTON SCATTERING APPROXIMATION 

In order to reduce (significantly) the computing time, we use an approximate Compton Scattering Kernel (CSK) that represents 
the physically important features of the exact CSK. We make the approximation that the scattering is isotropic in the plasma 
frame, and write the differential cross section as 

""' (v,n->i/,n') = ^-a c (v,T)f d ( v ,Ty). (Bi) 



dvdfl An 



Here, a c is the total cross section given in eq. d56h and fd is the spectral redistribution function of the photons. We require 
scattering to conserve photon number and require fd to satisfy 

/ f d (f>,f,i>')di>' = l, (B2) 

Jo 

and 

/■OO 

fd (v,f ,£>') vdv=va(v,f), (B3) 



10 

where vq is the average frequency of scattered photons. The approximations used for vq and fd are given below. We use different 
approximations for low (NR) temperatures and for high (relativistic) temperatures, with a transition temperature T„, = 0.25. We 
use a smooth interpolation between the two temperature regimes (over a ~ 10% interval in T). 

Low T (f < 0.25) 

Average energy shift — We chose vq to produce the correct average energy shift for v <C AT and for v ^> AT, and no energy shift 
for NR Compton equilibrium, v = AT. We use 



^o 

— = mm 

v 



47\4r+l)-i>(i>+l) \ AT 



(l+a u i>y 



(B4) 



for v < AT, and 

VIZ 



* i+ Mm 



(B5) 



for v > AT. a u {f) is determined by requiring that for a Wien spectrum, the energy gain of photons with energy less than AT, 



AT 





be equal to the energy loss of higher energy photons, 



Pgain{f,a) oc / dvv 2 e vlT <j c (v, f) (y (v, f ,a v ) - z>) , (B6) 



PlossChot fdi>D 2 e-^ f a c (i),i)(MvJ)-v)- (B7) 

if 
We use a 4-th order polynomial for a(log T), 

a v (f) = -0.003763 log(f) 4 -0.02311og(f ) 3 -0.019221og(f ) 2 -0.1291og(f ) + 3. 139, (B8) 

which is accurate to better than a percent, and set a v (T < 0.01) = a v (T = 0.01). 

Photon redistribution — We choose a photon re-distribution function that follows the shape of a thermal spectrum with a target 
temperature T lar = i>o(v, T)/A, 

fd{vJii>')=Av'\- 'l f "% (B9) 

where 
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High T(T> 0.25) 

High v - Klein Nishina corrections — For v > 1 /(AT), in the Klein-Nishina regime, we choose 

v a e-°'l f o(v' ,f) 



fd(PJ,#Y- 



a d (T)T 3 a T 



(Bll) 



The value of a d is chosen so that the integral over fd(v, f, v')dv' is 1 . This form ensures that a Wien spectrum with a relativistic 
temperature is unchanged by scattering of electrons with the same temperature. We use a 4th order polynomial approxiomation, 



a d (f) = -0.00461 llog(f) 4 + 0.0071971og(f) 3 + 0.09079 log(f) 2 -0.31661og(f) + 0.3146 



(B12) 



and set a(f) = a(5)(f/5) lJ for f > 5. This approximation is accurate to better than a percent for temperatures below m e c 2 , and 
to better than 10% everywhere. 

Low v - Inverse Compton — In order that a power law spectrum of the form I v oc v 1 retains its form after scattering, and in order 
to reproduce the ultra relativistic limit of the energy boost, 16T 2 , we choose 



f d (v,T,v')<x Vj7 



\n„ 



Q(v'-v). 



(B13) 



We use a cutoff at 8 1 to avoid overproducing photons at high frequencies, and normalize accordingly. 

C. PLASMA SPEED OF SOUND 

Below is a short derivation of a general formula for the speed of sound in a plasma of electrons, protons and <? + e" pairs, 
neglecting the thermal pressure of the protons (valid for T <C m p c 2 ). Here n\ is number density of leptons (electrons+positrons), 
and n p = «//(2x++ 1). For covenience we use p = p/m e c 2 , e = e/m e c 2 . Let s be the entropy per lepton. Rewriting Eqs. d38l and 
i we have 



p = niT, 



e = m[ l+f(T) + 



2x++l 



where we define 



where f(f) was defined in Eq. 
implies 



We therefore have 



and 



f(T)=-f(T)T, 



(CI) 

(C2) 

(C3) 



The derivative with respect to f at constant s is denoted by '. Constant entropy per lepton 



n, 






p =nif' + m 

^ + ( 1+ /+^)f. 



The speed of sound is finally given by 



P ss = c ss /c =\/'— = VT. 
e' 



1 + 1//' 



'-+=^ + f' 



(C4) 
(C5) 

(C6) 

(C7) 

(C8) 



This equation can be easily verified to obey the asymptotic NR and ultra relativistic limits. 

D. FRAME TRANSFORMATIONS 
Below are some useful transformation rules relating the values of v, I, r\ and \ measured in the shock and plasma rest frames: 



fJ'sh = 



M+/3 
1 + /V 



M = 



1 - PUsh 



(Dl) 
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I(V, v) ( v_ 

VshifJ-sh^sh) \Vsh) 



Xsh(Psh,V s h) 



V s h 



(D2) 
(D3) 

(D4) 

(D5) 
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